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Abstract 



In this thesis we study quantum transport properties of finite periodic quasi-one-dimensional 
waveguides whose classical dynamics is diffusive. We focus in the semiclassical limit which 
enable us to employ a Random Matrix Theory (RMT) model to describe the system. The 
requirement of normal classical diffusive dynamics restricts the configuration of the unit cells 
to have finite horizon and the appropriate random matrix ensembles to be the Dyson circular 
ensembles. The system we consider is a scattering configuration, composed of a finite periodic 
chain of L identical (classically chaotic and finite-horizon) unit cells, which is connected to 
semi-infinite plane leads at its extremes. Particles inside the cavity are free and only interact 
with the boundaries through elastic collisions; this means waves are described by the Helmholtz 
equation with Dirichlet boundary conditions on the waveguide walls. Therefore, there is no 
disorder in the system and all scattering is due to the geometry of the chain which is fixed. 
The equivalent to the disorder ensemble is an energy ensemble, defined over a classically small 
range but many mean level spacings wide. The number of propagative channels in the leads is 
iV and the semiclassical limit is achieved as N — > oo. An important quantity for the transport 
properties of periodic chains is the number of propagating Bloch modes Nb of the associated 
unfolded infinite periodic systems. It has been previously conjectured that for strongly diffusive 
systems in the semiclassical limit (N B ) ~ y/ND, where D is the classical diffusion constant. 
We have checked numerically this result in a realistic cosine-shaped waveguide with excellent 
agreement. Then, by means of the Machta-Zwanzig approximation for D we obtained the closed 
form expression (N B ) = y/N/ir, which agrees perfectly with the circular ensembles. On the 
other hand, we have studied the (adimensional) Landauer conductance g as a function of L and 
N in the cosine-shaped waveguide and by means of our RMT periodic chain model. We have 
found that (g(L)) exhibit two regimes. First, for chains of length L < \/N the dynamics is 
diffusive just like in the disordered wire in the metallic regime, where the typic ohmic scaling is 
observed with (g(L)) = N/(L + 1). In this regime, the conductance distribution is a Gaussian 
with small variance (such that (g^ 1 ) ~ 1/ (g)) but which grows linearly with L. Then, in longer 
systems with L 3> \/N, the periodic nature becomes relevant and the conductance reaches a 
constant asymptotic value (g(L — > oo)} ~ (Nb)- In this case, the conductance distribution loses 
its Gaussian shape becoming a multimodal distribution due to the discrete integer values N B can 
take. The variance approaches a constant value ~ y/N as L — > oo. Comparing the conductance 
using the unitary and orthogonal circular ensembles we observed that a weak localization effect 
is present in the two regimes. Finally, we study the non-propagating part of the conductance in 
the Bloch-ballistic regime, which is dominated by the mode with largest decay length £ which 
goes to zero as g np = 4e~ 2L//£ as L — > oo. Using our RMT model we obtained that under 
appropriate scaling the pdf P(£) converge, as N — > oo, to a limit distribution with an algebraic 
tail P(£) ~ £~ 3 for £ — > oo; this allowed us to conjecture the decay (g np ) ~ L~ 2 which was 
observed in the cosine waveguide. 



II 



Resumen 



En esta tesis estudiamos propiedades de transporte cuantico en gui'as de onda finitas per- 
iodicas quasi-unidimensionales, cuya dinamica clasica asociada es difusiva. Nos enfocamos en 
el lfmite semiclasico el cual nos permite emplear un modelo de Teoria de Matrices Aleatorias 
(TMA) para describir el sistema. El requisito de difusion normal de la dinamica clasica restringe 
la configuracion de la celda unitaria a tener horizonte finite, y significa que los ensembles apropi- 
ados de TMA son los ensembles circulares de Dyson. El sistema que consideramos corresponde 
a una configuracion de scattering, compuesto de una cadena finita de L celdas unitarias (clasi- 
camente caoticas y con horizonte finito) la cual esta conectada a dos gufas planas semi-infinitas 
en sus extremos. Las partfculas dentro de esta cavidad son libres y solo interactuan con los 
bordes a traves de choques elasticos; esto significa que las ondas son descritas por una ecuacion 
de Helmholtz con condiciones de borde tipo Dirichlet en las paredes la gufa. Por lo tanto, no hay 
desorden en el sistema y el scattering es debido a la geometrfa de la cadena la cual es estatica. 
El analogo al ensemble de desorden es un ensemble de energfa, definido sobre un intervalo clasi- 
camente pequeno pero cuyo ancho es varias veces un espaciamiento de niveles promedio (mean 
level spacing). El numero de canales propagativos en las gufas planas es N y el lfmite semiclasico 
se alcanza cuando iV — > oo. Un numero importante para las propiedades de transporte en ca- 
denas periodicas es el numero de modos de Bloch N B del sistema extendido infinito asociado. 
Previamente, ha sido conjeturado que en sistemas fuertemente difusivos en el lfmite semiclasico 
(Nb) ~ V ND, donde D es la constante de difusion clasica. Hemos comprobado numericamente 
este resultado en una gufa de ondas con forma de coseno obteniendo excelente concordancia. 
Luego, mediante la aproximacion de Machta-Zwanzig para D obtuvimos la expresion analftica 
{N B ) = ^N/ir, la cual concuerda perfectamente con los ensembles circulares. Por otro lado, 
hemos estudiado la conductancia (adimensional) de Landauer g como funcion de L y N en la 
gufa-coseno y mediante nuestro modelo RMT para cadenas periodicas. Hemos encontrado que 
(g(L)) muestra dos regfmenes. Primero, para cadenas de largo L < y/N la dinamica es difusiva 
tal como en un cable desordenado en el regimen metalico, donde se observa el escalamiento 
ohmnico tfpico con (g(L)} = N/(L + 1). En este regimen, la distribucion de conductancias es 
Gaussiana con una varianza pequeha (tal que (g^ 1 ) ~ I /(g)) pero que crece linealmente con 
L. Luego, para sistemas mas largos con L ^> y/N, su naturaleza periodica se hace relevante y 
la conductancia alcanza un valor asintotico constante (g(L — > oo)) ~ (N B ). En este caso, la 
distribucion de la conductancia pierde su forma Gaussiana convirtiendose en una distribucion 
multimodal debido a los valores discretos (enteros) que N B puede tomar. La varianza alcanza 
un valor constante ~ \/N cuando L — > oo. Comparando la conductancia para los ensembles 
circulares unitario y ortogonal, mostramos que un efecto de localizacion debil esta presente en 
ambos regfmenes. Finalmente, estudiamos la parte no-propagativa de la conductancia en el reg- 
imen Bloch-balfstico, la cual esta dominada por el modo con la longitud de decaimiento mayor £ 
que va a cero como g np = Ae~ 2L ^ e cuando L — > oo. Usando nuestro modelo de TMA obtuvimos 
que bajo un escalamiento apropiado la pdf P(£) converge, cuando N — > oo, a una distribucion 
lfmite con cola algebraica P(£) ~ £~ 3 para £ — > oo; esto nos permitio conjeturar el decaimiento 
(fl'np) ~ L~ 2 , el cual fue observado en nuestra gufa de ondas coseno. 
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Chapter 1 

Introduction 



The study of wave transport phenomena is quite an old subject dating back at least to 
the days of Newton, d'Alembert, Euler, Laplace and Bernoulli among others, when the wave 
equation itself was originally derived. Its first applications in the 1700's were to mechanical 
(elastic) waves such as vibrating strings and membranes, sound, water waves and later to light. 
In 1862, Maxwell showed that electromagnetic radiation follows the same equation. In 1927, 
Schrodinger proposed his equation to describe non-relativistic particles, which is a generalized 
wave equation with quadratic dispersion relation, where the field takes complex values and the 
potential energy plays a role similar to the internal index of refraction in dielectric media. The 
ubiquity of the wave equation in so many different physical contexts can be explained by the fact 
that it can be derived from a Lagrangian variational principle, where the system consists of a 
continuum of coupled harmonic oscillators. Since generic systems perturbed near an equilibrium 
point can be described to first order by a harmonic oscillator, the wave equation must hold at 
least to first order in the field amplitude. 

With the introduction of electromagnetic resonant cavities, waveguides and optical fibers, 
the properties of spatially constricted waves gained attention. For instance, microwave res- 
onators are closed metallic structures that confine electromagnetic waves with wavelength in 
the microwave part of the spectrum, this is from around one meter to one millimeter. Opti- 
cal cavities, where light between the infrared and ultraviolet part of the spectrum is confined 
with high reflective mirrors, are key components of lasers. Elastic resonators and waveguides 
have also been studied. In the realm of condensed matter, technological advances in the last 
decades have opened the possibility to manufacture devices at the mesoscopic scale, usually 
composed of metallic or semiconductors two-dimensional layers at the nanometer scale. For 
instance, the study of electric conduction in thin and highly disordered metallic wires in this 
scale brought about new physical phenomena such as weak localization and universal conduc- 
tance fluctuacions |Bee97] . which are a consequence of the coherent propagation of waves in a 
disordered medium. Other notable examples are GaAs-AlGaAs heterojunctions, where a thin 
conducting layer is formed at the interface between GaAs and AlGaAs[Dat95j, within which 
the electron's dynamic is well approximated by an effective mass description forming a two- 
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dimensional degenerate electron gas (2DEG). What distinguishes a 2DEG in GaAs is its low 
scattering rate due to impurities, meaning that at temperatures low enough to suppress phonon 
scattering the electron wavefunction coherence can be maintained at scales larger than the de- 
vice size and therefore its dynamic depends only on scattering due to boundary conditions in 
the confining two-dimensional layer, which can be designed at will in the laboratory. 

Spatially confined waves can be decomposed in a finite number of stationary modes which are 
associated to the solution of an eigenvalues problem where the operator is the time-independent 
wave equation. The shape of the cavity where the field is confined plays a major role in the 
statistical properties of this spectrum as well as in the solvability of the wave equation. In the 
small wavelength limit, waves are well described by geometric optics, this is by the propagation 
of rays which follow the classical trajectory of free particles (with possible position dependent 
speed if the potential or refraction index is not constant). Thus, its behavior depends crucially 
on whether the underlying classical mechanics is integrable or chaotic. A dynamical system 
with n degrees of freedom is integrable when it has n constants of motion; in this case the 
action-angle variables exist and in principle an analytical solution in term of quadratures can 
be found. When an integrable system is separable[Gut90j, the associated wave equation is 
also separable. Only systems with a high degree of symmetry are integrable; on the other 
hand, generic systems are non-integrable. Such systems are called chaotic and are characterized 
by the exponential divergence of nearby initial conditions in phase space. Therefore, closed 
analytical solutions do not exist for the classical dynamics in chaotic systems and the wave 
equation is not separable. Besides non-separability, one of the signatures of chaos in the wave 
equation is found on its spectrum. In chaotic systems, the statistical properties of the energy 
spectrum are universal after appropriate system-dependent scalings. This was first proposed 
by Wigner when studying the energy levels of heavy nuclei and later Bohigas, Giannoni, and 
Schmit[BGS84j conjectured the universality for general chaotic systems. Note, however, that 
the quantum (wave) dynamics in a classically chaotic system does not displays a property like 
exponential divergence of similar initial fields, because the wave equation (and Schrodinger 
equation) is linear in the field. Nevertheless, if one considers the evolution of an initial field 
under the action of two slightly perturbed propagators they diverge exponentially in the sense 
of expectation values [PZ02| . It is also possible to consider other characterizations of quantum 
chaos, for instance, in terms of the wave equation eigenstates[RS94, BSS98J. However, probably 
the most general and profound characterization of quantum chaos is the non-separability of the 
wave equation, which is shared with the Hamilton-Jacobi equation[Gut90j, allowing a direct 
connection between the quantum and classical domains. 
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1.1 Quantum dots 

One of the most common systems where quantum chaos has been studied, both in the 
theoretical and experimental literature, is the quantum dot. Quantum dots are confining cavities 
on the nanometer scale, usually carved in semiconductive material, were waves are only scattered 
by the geometric features of its enclosing walls; electrons in the cavity form a 2DEG. The 
underlying classical dynamics corresponds to a billiard where free non-interacting particles travel 
with constant velocity between elastic bounces against the quantum dot boundaries. Typically, 
the billiard dynamics is chaotic, so an ensemble of initial conditions will cover the full phase space 
volume after the ergodic time r e . In order to perform transport measurements, the quantum 
dot is attached to two electron reservoirs called contacts. If this coupling is weak so that the 
residence time in the quantum dot r r ^> r e , the spectral statistical properties are insensitive 
to microscopic details (in particular, to the geometry of the cavity boundaries). By applying a 
small potential difference to the contacts it is possible to measure the quantum dot conductance 
G, which at low temperatures and in the absence of impurities can be calculated using the 
Landauer formula |Lan57t IBILP85] , 



where e is the electron charge, h is Plank's constant and t is the transmission matrix. Landauer 
conductance is not related to an intensive resistance as is usually defined for macroscopic bulk 
conductors because it does not arise from phonon or impurity scattering; instead, one must 
think of a mesoscopic conductor (in particular quantum dots) as a complete phase-coherent 
unit, where resistance and transmission are a consequence of wave scattering and interference 
effects. All measurable scattering and transport properties of a quantum dot are encoded in the 
scattering matrix. 

Quantum dots display features similar to disordered wires shorter than its Anderson local- 
ization length, in particular weak localization and universal conductance fluctuations. Weak 
localization, which is the enhancement of reflection probability in time-reversal-symmetric sys- 
tems, was first considered in disordered wires by Abrahams et q/. |AALR79| and followed by ex- 
tensive theoretical and experimental work[Ber84j. Then, the universal conductance fluctuations 
(reproducible conductance fluctuations of order e 2 /h as a function of fermi energy or magnetic 
field) were theoretically predicted in disordered wires by Altshuler[BL85j and Stone [Sto85j, and 
observed experimentally for the first time by Webb and Washburn |WWUL85j in metallic rings. 
Several years later, Baranger, Jalabert and Stone|BJS93j developed a semiclassical theory to pre- 
dict these effects in quantum dots which were also observed experimentally [MR W + 92| lBIA + 94] . 
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1.2 Semiclassical theory 

The first efforts to understand the effects of an underlying chaotic classical dynamics in the 
quantum level were focused on closed systems, in particular to the study of the energy spectrum's 
statistical properties. In this regard, Gutzwiller put forward his celebrated trace formula[ Gut71l 
IGut90j . a semiclassical expression for the density of energy states in terms of a sum over classical 
periodic orbits. Subsequently, Gutzwiller theory was generalized by Eckhardt et al. |EFMW92] to 
express arbitrary quantum observables matrix elements in terms of classical periodic orbits. On 
the other hand, the work of Berry and Tabor [BT77J and Bohigas, Giannoni and Schmit[BGS84] 
focused on the universal character that many empirical and numerical evidence attributed to the 
spectrum of chaotic systems. They pointed out the fact that classically integrable and chaotic 
systems had quite different level statistics in the semiclassical limit; in particular they proposed 
that the nearest-neighbor level spacing distribution of integrable systems follows a Poisson law 
whereas chaotic systems show a higher degree of repulsion with a distribution similar to the one 
implied from ensembles of random hamiltonian with Gaussian elements. This is the foundation 
of modern Wigner-Dyson random matrix theory for Hamiltonians in chaotic systems and one 
of the better understood manifestations of quantum chaos. Recently, starting from the trace 
formula, it was proved that the Wigner-Dyson random matrix theory correctly describes the 
spectral statistics of quantum chaotic systems in the semiclassical limit |MHB + 041 |HMA + Q7] , 
i.e. the Bohigas, Giannoni, Schmit conjecture is fulfilled. 

In closed chaotic systems there is only one dimensionless parameter relevant for its asymp- 
totic semiclassical properties which can be taken as the wavelength in units of the system typical 
size, X/L. In the limit X/L — > the universal semiclassical regime is reached and energy fluc- 
tuations are described by random matrix theory. Note that this is formally equivalent to the 
limit ?z. — >■ which is a common way to take the semiclassical limit. On the other hand, the 
semiclassical treatment of open chaotic systems is more complex due to the existence of an addi- 
tional scale given by the mean residence time r r , which in case of closed systems is infinite. The 
exponential divergence of nearby initial conditions in chaotic systems implies that the quantum 
evolution of an initial wave packet follows closely the classical dynamics up to the Ehrenfest 
time, tE ~ A -1 In S/h [Ber79j, where A is the sum of positive Lyapunov exponents and S is 
a characteristic classical action such as that of the shortest periodic orbit. Thus, a universal 
statistical description of open chaotic systems in terms of random matrix theory only holds 
for tE <C T r and a system-dependent regime is observed for r r <C Note that for quantum 
dots, the requirement of satisfying the limit r e <C r r means that they fall in the random matrix 
universal regime for fixed h. A semiclassical theory for quantum scattering in open chaotic sys- 
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terns was first worked out by Baranger, Jalabert and Stone[BJS93j in terms of a semiclassical 
propagator, similar to which had been used before by Gutzwiller for the trace formula. They 
wrote the t matrix elements t nm for a two-terminal quantum dot as a sum over trajectories 
connecting mode m on the left entry to mode n on the right. This theory correctly predicts the 
weak localization and universal conductance fluctuations observed in quantum dots. A random 
matrix theory for transport was then introduced by Jalabert, Pichard and Beenakker[JPB94j 
using the Dyson circular ensembles to model the quantum dot scattering matrix. 



1.3 Diffusion in extended chaotic systems 

It is well known that in extended classical chaotic systems, under general assumptions, the 
dynamics relax by a diffusion process fGas98| . This means that an initially confined ensemble of 
initial conditions will spread as (ccf) — > Dt for t — > oo, where x t is the position vector at time 
t, D is the diffusion coefficient and the average (•) is taken over the equilibrium measure. One 
of the earliest theories of transport is due to Drude, who tried to explain electron conduction in 
metals from a purely classical perspective. An equivalent, more contemporary formalism, is the 
Green-Kubo linear response theory |Dor99) . Drude's model treats electrons as non-interacting 
(but electrically charged) particles traveling in a lattice of infinite mass spheres representing 
the ions; the only interaction are elastic collisions between particles and spheres. When a small 
electric field E is applied to the metal, the electric current is given by Ohm's law J = aE, 
where the conductivity can be written as a function of the electron density of states at the Fermi 
energy, p(E F ), by 

a = e 2 p(E F )D. (1.2) 
This is known as the Einstein relation for degenerate conductors. The quantum property of 



electrons being distributed according to a Fermi-Dirac distribution is implicit in equation (1.2); 



all other quantum effects are neglected. In this sense, (1.2) is a semiclassical result. Regardless 
of its simplicity, Drude's model provides a good approximation for the electric conductivity, 
Hall effect and thermal conductivity in metals at room temperature. 

The quantity which can be really measured experimentally is not the conductivity but the 
conductance G = I /V, with / the electric current and V the applied voltage. In a macroscopic 
homogenous two-dimensional conductor, Ohm's law implies that 

G=^, (1.3) 
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where W and L are the conductor width and length, respectively. Obviously this scaling must 
break down when we go to the mesoscopic scale and quantum phase coherence starts playing 
a role. In particular, the conductivity loses its meaning for systems smaller than the phase- 



relaxation length, and in this case, the conductance is given by the Landauer formula (1.1) in 
terms of transmission probability. One of the most striking novel features observed in meso- 
scopic conductors is Anderson localization in disordered wires, which consists in the exponential 
suppression of transmission (and hence of conductance) in a disordered wire longer than the 
localization length £. Notably, however, for wires of length in the range £ < L < £, with I the 
mean free path, the conductance displays ohmic scaling with G ~ L" 1 ; this is the metallic or 
quantum diffusive regime, in which weak localization and universal conductance fluctuations are 
observed. Edwards and Thouless noticed that conductivity in disordered systems is closely re- 
lated to the sensitivity of its eigenvalues on an external perturbation |ET72] and argued that the 
metal-insulator transition could be described by a single parameter, namely the ratio between 
diffusion energy E c = hD/L 2 (which is the energy associated to the time it takes a classical 
particle to escape the systems by diffusion) and the mean level spacing AE. In fact, the quan- 
tity g c = E C /AE has been shown to be the average conductance in disordered systems in the 
metallic regime |SA93bJ. For short systems in the semiclassical limit, the conductance starts 
being g c ^> 1 and decays as until g c ~ 1 when Anderson localization kicks in. 

Localization is a purely quantum (wave) effect, caused by the destructive interference of 
multi-scattered waves inside disordered media. In the opposite extreme is the case of periodic 
media, where there is other well known quantum result, namely Bloch theorem, which states 
that in infinite periodic systems the energy spectrum is organized in bands and eigenstates 
are ballistic, i.e. they travel through the system with constant velocity. This means that 
transport in extended periodic systems is ballistic even if the underlying classical dynamics 
is chaotic. The signatures of classical diffusion in band spectra where studied by Dittrich et 
a/.[DMSS97, DMSS98J using the semiclassical trace formula and by Simons and Altshuler[SA93bJ 
in a more general context using the non-linear a model. Simons and Altshuler showed that the 
energy levels of chaotic and disordered systems subject to an external perturbation (which for 
periodic systems is identified with an Aharonov-Bohm flux) are universal up to two system 
dependent parameters. Later, Dittrich et al. studied the two-point correlation function in finite 
periodic systems using semiclassical methods and found agreement with the results of Simons 
and Altshuler in the large system limit. 

In this thesis we are interested in the semiclassical transport properties of periodic systems 
possessing an underling chaotic or diffusive classical dynamics. We focus in two related quanti- 
ties, the number of propagating Bloch modes Nb and the (adimensional) Landauer conductance 
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g. We show that in periodic chains, (g) ~ (Nb) in the long-chain limit. An analytical expres- 
sion for (Nb) was calculated by Faure[Fau02j in the semiclassical limit. We test this result in a 
physically realistic cosine waveguide model, show its connection to random matrix theory and 
obtain an universal expression for (Nb) in a infinite chain of quantum dots with and without 
time-reversal symmetry, revealing a weak-localization-like effect in this quantity. On the other 
hand, we study the conductance properties of diffusive periodic systems employing a cosine 
shaped waveguide and a random matrix quantum dot chain model. We find that an ohmic 
regime is observed in such systems, but over a shorter length range than in disordered systems. 
The conductance can be separated in a Bloch-ballistic term, which remains non-null in the long 
system limit, and a decaying part related to closed modes. We show that the non-propagating 
part of the conductance decays with a power-law as L~ 2 . Finally, we consider the conductance 
distribution and fluctuations as a function of length and how the diffusive to Bloch-ballistic 
transition is signaled in them. 

1.4 Outline 

This thesis is organized as follows. In chapter [2] we review the scattering formalism for two- 
dimensional waveguides and Bloch theorem in such systems. By means of Oseledets theorem we 
show that the (adimensional) Landauer conductance is bounded above by Nb- Additionally, we 
show formally that the conductance is quasi-periodic as a function of length in the Bloch-ballistic 
regime. In chaper |3j we present relevant aspects of chaotic classical dynamics in billiards; 
in particular, we write an explicit result for the diffusion coefficient in the Machta-Zwanzig 
approximation for arbitrarily shaped quantum dots. The billiard model we use in the next 
chapters, namely the cosine billiard, is studied numerically. Chapter [4] is a digression on random 
matrix theory and its applications to study chaotic systems. We define Dyson circular ensembles 
and present the quantum dot periodic chain model. In chapter [5] we review Faure's result for 
(Nb) and show numerical calculations done with the cosine waveguide. A connection with the 
universal parametric correlation function of Simons and Altshuler is presented. We also derive a 
parameterless expression for (Nb) in a quantum dots chain, and using this model we characterize 
some of the features of the Nb probability distribution. The main results regarding transport 
properties are in chapter [6| where we study the conductance in the periodic cosine waveguide 
and random matrix models. First, we focus on the Bloch-ballistic regime, and then on the 
ohmic regime in periodic diffusive systems. We also show the existence of a weak localization 
correction in both regimes and study the conductance fluctuations. Finally, in chapter [7| we 
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conclude and present future work prospects. 



8 



Chapter 2 

Waveguides 



In this chapter we give an account of the quantum and classical dynamics of a particle 
inside a waveguide. The most general definition of a waveguide is a structure able to confine 
and guide waves. They can be found in a number of contexts such as electromagnetism - with 
different application on different parts of the spectrum- and acoustics. More recently, meso- 
scopic waveguides have been constructed where electrons are constricted in coherent metallic 
(or semiconductor) rings and wires with transversal section a few effective electron wavelengths 
wide. Depending on the underlying physical system, the wave equation and boundary condi- 
tions defining the waveguide mathematically can vary. We consider a two dimensional system 
where the wave dynamics is governed by the Helmholtz equation 

(V 2 + A; 2 )^ = 0, (2.1) 

for the field \1/ = ty(x,y), which is equivalent to the Schrodinger equation for a free particle 
with energy E = h 2 k 2 /2m. The boundary conditions are Dirichlet f = on the waveguide 
borders. This means that waves propagate coherently in the guide and are only scattered by 
the geometry of the boundaries. The associated classical dynamics to this Schrodinger equation 
is given by free non-interacting point particles colliding elastically with the waveguide walls; 
such systems are usually called billiards and are discussed in more detail in chapter |3j In this 
work we are interested in waveguides whose classical dynamics is strongly chaotic, so an initially 
bounded ensemble of particles will spread diffusively in the system such that the average square 
displacement (Ax 2 ) t ~ t with time t — > oo. For the wave equation (and Hamilton- Jacobi 
equation), classical chaotic dynamics translates into non-separability and hence, in general, the 
lack of analytic solution. Note that it is also possible to describe charged particles subject to a 



magnetic field with equation (2.1) imposing minimal coupling, V — > V — ieA/h where A is the 
vector potential and e the particle charge. Since the wave equation we study is the Helmholtz 
equation, the results we obtain are also valid for classical waves. 



In section 2.1.1 we review how to describe the scattering properties of a waveguide by means 



of the scattering and transfer matrices, and how to obtain its Landauer conductance from the 
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transfer matrix spectrum using the polar decomposition. In addition, we define the Bloch 

we 



spectrum of periodic waveguides using the scattering approach. Finally, in section 2.2.1 



obtain an expression for the conductance of a long periodic waveguide using Oseledets theorem. 



2.1 Scattering and transfer matrices 

2.1.1 Definition and properties 

A natural way to describe the wavefunction *&(x, y) in the waveguide is to project it on the 
local transverse basis, this is, writing the wavefunction as 



(2.2) 



71=1 



where g n (x,y) are the local transverse modes which satisfy the boundary conditions on each x, 
and c+(x) (c~(x)) is the right-going (left-going) longitudinal mode. In the particular case of a 
hard- wall waveguide, 



9n(x,y) 



h(x) 



sin hit 



V ~ h(x) 
h 2 (x) - hi(x) 



n 



oo, 



(2.3) 



where h\(x) < h 2 (x) are the wall heights as a function of the longitudinal coordinate x. This set 
of functions satisfies the null boundary conditions g n (x,hi(x)) = g n (x,h,2(x)) = everywhere 



in the guide. The longitudinal modes are obtained by inserting (2.2) in (2.1), which transforms 



the original partial differential equation into a system of coupled ordinary differential equations 
which can be efficiently solved numerically. 



In a plane lead, h\(x) and h 2 (x) are constant, so (2.1) is separable since g n (x,y) = g n (y) is 
independent of x. In this region, the longitudinal modes are given by 



e n( x ) 



(2.4) 



where k^ = k 2 — {n-n/ho) 2 is the longitudinal wavenumber, h Q = h 2 — hi is the lead width and 
the normalization is to impose unit flux. There are 



N 



h n k 



71 



(2.5) 
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A 1 


S matrix 


A r 
► 


► 




* B> 




B r 



Figure 2.1: Schematic diagram of the wavefunction decomposition in the plane leads. The 
gray area represents the scattering zone with scattering matrix S. The far field is completelly 
represented by the vectors A 1 , B l , A r and B r . 



propagating modes because for n > N the longitudinal momentum hk n is imaginary, implying 
null energy flux. These are called evanescent modes and decay exponentially with x. 



We will consider a scattering geometry [see figure 2.1 , i.e. a waveguide composed of an 
interaction or scattering region connected to two semi-infinite plane leads at its extremes. The 
far field wavefunction in the leads can be described with a 2N dimensional complex vector 
composed of coefficients A n and B n for n — 1, . . . , N and we can write the field as 

N 

*(x,y) = (^e+(x) + B n e~{x)) g n {x,y). (2.6) 

71=1 

Let A T (B v ) be the N dimensional complex vector of right-going (left-going) amplitudes in the 
right lead and A 1 (B l ) the same on the left lead. We denote the incoming and outgoing (or 
incident and scattered) fields in vector notation as 



V'in = „ T and V out = . (2.7) 





The scattering matrix S is defined as the linear transformation that maps incoming to outgoing 
fields, 

</w = Sip- m . (2.8) 

If t and r (tf and r') are the left-to-right (right-to-left) transmission and reflexion matrices, we 
have 

rA + t'B 1 = B l 



thus, we can write S as 



(2-9) 

tA l + r'B v = A\ 



(2.10) 
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which is a 2N x 2N dimensional matrix. The conservation of probability flux (or energy flux 
for classical waves) tp\ n ip^ n = ip OVLi ipout implies that S is a unitary matrix, i.e. 

&S=1, (2.11) 

with the hermitian conjugate of A. On the other hand, if time reversal symmetry is not 
broken by a magnetic field the scattering matrix is symmetric, 

S = S T , (2.12) 

with A T the transpose of A. 

The archetypal waveguide configuration we consider in this thesis will be a periodic chain 
of scatterers. Given the unit cell scattering matrix Si and assuming the coupling of evanescent 
modes between cells in the chain is negligible, it is possible to obtain the (n + l)-cells chain S 
matrix by the usual Feynman paths sum rules[Dat95j, 

r n+ i = r n + t' n r \(1 - r^r^Hn 
t n +l = - r'^xYHn 

(2.13) 

^ + l = ^+*l^n( 1 -nO~V 1 

where the subscripts denotes the number of cells in the chain. If the coupling of evanescent 
modes between cells cannot be neglected, transmission and reflexion matrices are in general not 
invertible and other methods to relate unit cell and chain scattering matrices is necessary [see 
appendix [A], 

We have seen that the S matrix links incident and scattered waves. A complementary 
picture is given by the transfer matrix which connects the wavefunction on the right and left 
leads, 

bright = ( ^ r j and Vieft = ( ^ j • (2.14) 
Then the transfer matrix M is defined by 

bright = MVtat- (2.15) 



12 



52.1 Scattering and transfer matrices 



Waveguides 



In order to obtain an explicit form for M we recast (2.9) to the form 



Alright = M r V>l e ft 



(2.16) 



with 




M,. = I I and Mi = | I .', j . (2.17) 

Hence, the transfer matrix can be written in terms of transmission and reflexion matrices as 




M = M[ l M r . 



(2.18) 



From the unitarity of S (2.11) follows that the M matrix satisfies 



MSM f = £, 



(2.19) 



(2.20) 



where 

ljv On 
On —In 

with ljv and 0^ the N x N identity and zero matrices. Thus, M is a pseudo-unitary matrix^ 
with metric S whose induced norm ip^Y^ip can be interpreted as the probability flux of a field 



ip through the system. In addition, the scattering matrix symmetry (2.12) in time-reversal 
symmetric systems implies[BDR10j a further symmetry of M given by 



KMK = M*, 



where 



K 




(2.21) 



(2.22) 



The assumption of negligible coupling between adjacent cells required for (2.13) implies that 



it is possible to truncate the unit cell reflexion and transmission matrices keeping only open 
channels without loosing significant accuracy. When this is the case we can relate the unit cell 
transfer matrix M\ to the length L periodic chain transfer matrix M by 



M = (Mi, 



(2.23) 



The pseudo-unitary group defined by a metric with a signature equivalent to S is usually denoted U(N,N). 
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This is formally equivalent to the S matrix composition equations (2.13) but numerically un- 



stable since Mi contains unbounded elements. In appendix [X] we show an alternative method 
to obtain the length L chain S matrix using the unit cell M matrix spectrum. The numerical 
integration procedure we employ to obtain the S matrix is due to V. Pagneux and is discussed 
in detail in |PaglO|. 

2.1.2 Bloch spectrum of periodic waveguides 

Lets considered the unfolded periodic waveguide, i.e. the periodic scatterers chain of infinite 
length. As is well known, Bloch theorem states that in a system invariant to (discrete) transla- 
tions in the x direction Schrodinger equation solutions take the form (j) nA {x,y) = e iqx u n ^ q (x, y), 
with u n , q (x,y) a periodic function of x with a period equal to the translation length (or equiv- 
alently the potential periodicity or unit cell length). Moreover, the energy levels of the system 
form continuos bands parametrized by the quasi-momentum q and the integer index n. One 
method to obtain the set of Bloch eigenfunctions and eigenvalues is solving the hamiltonian 
eigenvalue problem directly in the unit cell domain for arbitrary quasi-momentum q. This 
means imposing periodic boundary conditions such that 



J n,q 



[x + Z,y) = e 



iqZ _ 



J n,q 



>, y) for all y, 



(2.24) 



with Z = x\ — Xq the unit cell length, xq < x < X\ and xq and x\ the unit cell borders. This 
produces a discrete number of solutions E = E n (q) and <p n ^ q (x,y) labeled by n. On the other 



hand, we see from (2.15) that we can impose Bloch condition (2.24) for a fixed energy E using 



the unit cell transfer matrix M\ (calculated for that energy) in the eigenvalue problem 



Mxij) = Xi[). 



(2.25) 



The advantage of this scattering approach is that once we have computed the S matrix it only 
takes the solution of a matrix eigenvalue problem to obtain the Bloch basis. In addition, using 



(2.16) we can rewrite (2.25) as the generalized eigenvalue problem 



(2.26) 



which is numerically more robust since it avoids taking the inverse of M\. 

Solving (2.25) produces 2N eigenvalues Aj = e iqiZ and associated eigenstates ipi, with qi in 



general complex numbers. Eigenvectors ipi associated to eigenvalues A, with unit modulus (or 
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equivalently real quasi-momentum q{) correspond to propagating quasi-periodic Bloch-Floquet 
modes whereas states with eigenvalues such that |Aj| 7^ 1 are closed channels in the unfolded 
infinite periodic system and have null flux = ip^Hip. Propagating Bloch modes are right-going 
or left-going if > or < respectively and have a group velocity given by 



V n (0) 



Z dE„ 



h dO 1 



(2.27) 



with 9 = qZ. Note that the velocity v n (8) is proportional to the energy band E n (6) slope at the 
intersections E = E n (8). Hence, since the bands must be 27r-periodic, the number of right-going 
and left-going Bloch modes must be the same. 



From the transfer matrix pseudo-unitarity (2.19) follows|HP07j that if ip is a right eigenvector 
of (2.25) with associated eigenvalue A, then (5/0)* is a left eigenvector with eigenvalue 1/A*. 



This means closed channels Bloch quasi-momentum comes in pairs (q, q*), which correspond to 
states decaying in both chain directions. In addition, for time-reversal symmetric waveguides 



we have from (2.21) that if ip, A is an solution of (2.25) then Kip, A* is also a solution, so 



in this case the M matrix eigenvalues come in quadruples (A, A*, 1/A, 1/A*). This means that 
open channel Bloch quasi-momentum comes in pairs (q, —q), i.e. the bands are symmetric with 
respect to q = when time-reversal symmetry is preserved. 

In finite periodic waveguides, the Bloch basis is also important since transmission in a long 
chain will only be allowed for modes that are linear sums of the elements of this subspace 
{ipi : I Aj| = 1}. Since M = {M\) L the length L chain transfer matrix spectrum is given by 



,Ll N 



{Xf }f = i, hence closed Bloch modes will decay exponentially in the chain. In section 2.2.2 
will see this in more detail. 



we 



2.1.3 Polar decomposition 



The unitarity of the S matrix (2.10) implies that the Hermitian matrices tv , t't'\ 1 — rr^ 



and 1 — r'r't have the same espectrum {Ti}f =l , with all elements real and bounded in the [0, 1] 
interval. The scattering matrix can be written in terms of this set of eigenvalues by means of 
the polar decomposition |PAPN 88j . 



U 
V 




(2.28) 
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where U, V, U' and V are unitary N x N matrices and T is a diagonal matrix with elements 
7i,T2, . . . ,T/v. The transfer matrix can also be written using this decomposition, 



M 



V 

v't 



7 r 1 - 1 



17' 
C7t 



(2.29) 



From (2.19), we have that the spectrum of MM^ consists of inverse pairs {Aj, Aj 1 }fL 1 , which 



are real and positive. From (2.28) and (2.29) one obtains 



[21 + MM f + (MM" 1 ") 



tWl-i 



1 I ttf 
4 1 t't't 



(2.30) 



which is a 2-to-l relationship between the spectrum of MM\ {Aj, A i and the spectrum 
of tt\ {Tj}f =1 , given by 

4 

T i = T , i = l,..., A. (2.31) 







Figure 2.2: Example of a periodic chain with four unit cells. We denote the unit cell length Z 
and the leads width h . 



2.2 Asymptotic conductance of a periodic waveguide 

2.2.1 Landauer conductance 

Landauer scattering theory of electron conduction [F L811 IBILP85) provides a complete 
description of transport in mesoscopic devices at low temperatures, voltages and negligible 
electron-electron interaction. The typical system consists of a phase- coherent region where 
scattering is elastic (the periodic chain in our case) connected by ideal leads to two electron 
reservoirs at zero temperature and energies E and E + 8E, respectively. The small potential 
difference SE E generates an non-equilibrium net current between the reservoirs, 

I=^-gSE, (2.32) 
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where e is the electron charge and g = tr[tt'] is the iV-channels transmission probability through 
the scattering zone. The conductance G is defined by I = G5E/e, hence 



G 



7T 



h 



9 



7T 



h 



tr[tt f 



(2.33) 



Note that if transmission were perfect tr[tv] = N, thus Landauer's conductance implies a non- 
null resistance G' 1 = 7ih/(Ne 2 ) even for reflectionless samples, known as contact resistance, 
which arise from non-elastic scattering when the discrete N channels propagating in the leads 
couple to the reservoirs continuous spectrum and thermalize at the equilibrium potential. This 
is a purely quantum effect since it vanishes when iV — > oo. 

Conductance seen as transmission is not only relevant for conduction in mesoscopic degen- 
erate electron gases but is an important quantity in other contexts when we are interested in 
understanding transport properties of waves, for example in electromagnetic and photonic de- 
vices, in molecular nanowires and in acoustic waveguides. In addition, is possible to define other 
transport quantities as a function of transmission eigenvalues Tj such as shot noise power and 
thermopower |B92 lKSG + 96l lGMB + 99] . For the sake of broadness and notation simplicity, in 
what follows we consider the adimensional conductance 



N 



(/ = tr[tt t ] = ^T i . 



(2.34) 



i=l 



2.2.2 Oseledets theorem for transmission eigenvalues 

The M matrix spectrum of a periodic chain is trivially obtained from the unit cell Bloch 



spectrum since as we have seen Aj(L) = Xf . However, the conductance (2.34) -and in general 



any other transport property dependent of the transmission eigenvalues Tj- is a function of 



MM* 1 eigenvalues as shown by (2.31). Thus, g depends on the chain length L in a non trivial 



way because of composition relations (2.13). Nevertheless, it is possible to obtain an asymptotic 



approximation for g from the unit cell transfer matrix and Oseledets theorem |Qse68j . Let {Aj} 
be the set of M\ eigenvalues, then Oseledets theorem implies that 



2N 
i=l 



lim (MfAff^sz 

L— >oo 










|A 2 | 



\ 





|A 2 at| 



(2.35) 
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this is, A i (L) 1 /( 2L) -+ | Ai| where {Aj(L), A i (L) _1 }^ 1 is the spectrum of MM f 
Hence, we have that 

Aj(-L) — ► a^e^H 

L— >-oo 



m[m[ x . 

(2.36) 



where a,(L) is a positive and (generically) bounded function of L. Then, using relation (2.31), 
we can decompose g in two terms, one with the sum of transmission modes related to the 2Nb 
propagating Bloch modes |A$| = 1 and another with the sum of modes related to evanescent 
Bloch states |Aj| 7^ 1, which have a decay length 



min{log|Ai|} 

|Ai]>l 



(2.37) 



determined by the slowest to decay non-propagating state. From (2.36) and (2.31) we deduce 



that the transmission eigenvalues Tj associated to Bloch modes is of order one, so for chains of 
length L>£, 

(2.38) 



^(L)<iY B + 4a m (L)- 1 e 



-1 p -2L/l 



where a m is the Oseledet function associated to I. The equality in (2.38) is non-generic 



and occurs if Mi is normal, in which case Tj(L) = 1 for all Bloch modes. In appendix [BJ a 
perturbative calculation for the Bloch transmission eigenvalues Tj(L) is performed in the limit 
Mi is near a normal matrix. For the generic case, g(L — > 00) ~ Nb and only once we had 
defined the appropriate ensemble averages for g and Nb in chapter [4] we will be in condition to 
establish an exact relation between this quantities. This is done in chapter [6j 



2.2.3 Conductance quasi-periodicity 



Not unexpectedly, (2.38) contrasts with quasi-one-dimensional disordered systems where 



waves localize leading to zero conductance in a long wire. In a periodic chain, if the associated 
unfolded periodic system has a non-null number of propagating Bloch modes, the conductance 
do not decay to zero and instead show a quasi-periodic behavior in the big L limit. We can see 
this as follows. From the transfer matrix definition (2.17) follows that t = (X^)^ 1 where 



M 




(Mi 



PD L P \ 



(2.39) 
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The matrix P = [i^]^ has the set of Mi right eigenvectors in its columns and P 1 = [uj 
the set of Mi left eigenvectors in its rows. Then, 



T12V 
i=l 



2N 



f 1 = E tf/W ( 2 - 4 °) 

2=1 

where A and r/j are iV dimensional vectors composed of and it; second halfes respectively. 



We separate the sum in (2.40) in two parts, one with the 2Nb unit modulus eigenvalues Aj and 
another with the rest, where in the long chain limit all terms with |Aj| < 1 can be neglected so 
iV — Nb elements survive. Then, 

t" 1 = A + UG L V, (2.41) 

where 

(2.42) 
(2.43) 
(2.44) 
(2.45) 

i.e. A is a iV x N matrix, U is a N x (TV — Nb) matrix with the vectors {A : |Aj| > 1} in 
its columns, V is a (N — Nb) x N matrix with the vectors {r/j : |Aj| > 1} in its rows, and 
G is a diagonal matrix with N — Nb elements {Aj}j ; |A 4 |>i- Note that A may not be invertible 
since, at best, it has rank 2Nb which can be less than N. To obtain t for L ^> £ using this 



A = 


Ei: 1*1=1 AiTW 


U = 


[A]i: |Ai|>l 


V = 


[ r ? i T ]i:|A l |>l 


G = 


diag({Ai} i: | Ai |>i), 



decomposition, we invert (2.41) using a generalized Woodbury Identity |Rie92] . which leads to 



t = A + - A + U { \C\ - C 2 Vj[A + + C 2 (G- L + V^A+uA C\, (2.46) 

where C\ = U± \ U±Ux) , C 2 = V± (v±V±^j an d A + is the Moore-Penrose generalized 
inverse of A. For any matrix Z, Z\\ = PaZ and Z± = (1 — Pa)Z with Pa the projector to the 
column space of A. Note that G~ L biggest element decay as eT 2L l l and U, V are independent 
of L. Finally, we have that for long chains g = tr[tt<] depends only on L quasi-periodically 
through the propagating Bloch modes in A + given by Mi eigendecomposition. 

The quasi-periodicity of g(L) can be understood as a consequence of the coupling between 
Bloch modes in the periodic chain and the plane leads propagating channels. The amplitude 
of this coupling, i.e. of the transmission between periodic chain and plain leads, oscillates as a 
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function of L because of the t matrix resonances due to interior reflection. 
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Chapter 3 

Classical dynamics in chaotic billiards 



The classical analog of the waveguide system presented in chapter [2] is an open billiard. 
An open (closed) billiard is a bidimensional system of non-interacting point particles enclosed 
in a infinite (finite) domain delimited by hard elastic walls. In open billiards, one is usually 
interested in the transport properties (relaxation) of an initially spatially confined ensemble 
of particles or in the scattering properties of an isolated interaction region. Of course, both 
problems are connected in the sense that transport properties can be thought to be a result of 
local scattering phenomena |LP67] . In chapter |2j we have defined the conductance of a periodic 
chain -which is a measure of transport probability- as a function of the scattering matrix of a 
single cell. This scattering approach to transport is also natural to experimentalists because, in 
the laboratory, finite pieces of materials are proved and then its properties extrapolated to the 
bulk. 

The classical dynamics of a billiard, depending on the shape of its boundaries, can be 
integrable, chaotic or mixed. Billiards have been extensively studied in the context of classical 
chaos -both in the physical and mathematical literature- because their simplicity has allowed 
several analytical results for particular billiards such as the Lorentz gas (Sinai billiard) and the 
Bunimovich stadium where, for instance, hyperbolicity, ergodicity and the central limit theorem 
have been proved. 

In this chapter we define a billiard in more detail and give a brief account of properties 



relevant for this thesis. In section 3.1 we define the billiard collision map and present the concepts 



of chaos and hyperbolicity. Then, in section 3.2 some statistical properties of deterministic 
chaotic systems are discussed. A good treatment of these topics can be found in |Gas98) . We 
describe the Machta-Zwanzig approximation in section 3^3 which will be used in chapter [4] to 
obtain an analytical expression for the number of propagating Bloch modes in a random matrix 
periodic chain model. Finally, the cosine billiard is defined; this is the system employed in our 
numerical calculations of the conductance properties discussed in chapter [6) 
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3.1 Billiard dynamics 

3.1.1 Billard flow and collision map 

A billiard is a hamiltonian system of non-interacting free particles bouncing between obsta- 



cles where they collide elastically [see figure 3.1 . We consider bidimensional billiards only; their 
hamiltonian is simply H = (p 2 . + p 2 )/2m so particles travel with constant speed. Without loss 
of generality we can set the particle speeds and mass to one and rescale later if necessary. Let 
V C M 2 be the position space where the particles can move. The billiard phase space is given 
by Ai = T> <g> S 1 where S 1 is the unit circle which represents all possible directions of motion. 
The continuos time evolution in the billiard consists of linear segments of free flights in T> and 
specular reflection at the boundaries &D. We can denote a phase space point in the billiard as 
X = (q,v) with q = (x, y) G T> and v G S 1 . Formally, the dynamical system is defined by the 
system of equations 

X = F(X), (3.1) 



where 



F(X) = £ • d x H 



(3.2) 



with X the symplectic fundamental matrix. Equation (3.1) induces a continuous time evolution 
or flow denoted by 

X t = $ ( X 0l (3.3) 
where X and X t are the phase space coordinates at time zero and t respectively. 




Figure 3.1: Generic billiard with convex obstacles (gray objects) and a particle trajectory. If 
all the dashed lines were part of the boundary 8Dq this would be a closed billiard. Otherwise, 
taking periodic boundary conditions on the dashed lines we could construct a periodic billiard 
with this unit cell. 



Let T(X) be the time of first collision of a phase space point X, i.e. the time X takes to 
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reach the first intersection with dT>. Given an initial condition X to = (qo, vq) we can iteratively 
define the billiard flow as 



^X tn = &~ T ( x ^x ti 

n„ ■ Vn)n„. ) 

(3.4) 



X tl = (q tl , v - 2(n ftl ■ v )n qti ) 

qtt = qo + v T(x to ) 

t 1 =t + T(X t0 ), 



where n q G S 1 is the vector normal to the billiard boundary at q G &D. The piecewise 
continuos nature of the billiard flow induces naturally the billiard map, which is the discrete 
time version of $ that iterates particles from collision to collision, i.e. is the Poincare map 
for the surface of section given by the billiard boundary. We define the billiard map for 
elements x = ( s ? sin 6) G dT> <g> [—1, 1] where sin# is the tangential momentum at the collision 
point r = r(s) with s the boundary curve natural parametrization; these are called Birkhoff 
coordinates. The map, 

Xn+l = 4>Xn , n . 

t n +l — tn + T{Xn)i 



is defined using (3.4) taking X to = (r n (s), cosdn rn +sin8t rn ) with t r the tangencial boundary 



vector h r ■ t r — 0. Note that there is ambiguity in the choice of t r direction; for instance, we 



can define it such that t r ■ v > for 9 > 0. See figure 3.2 

For open periodic billiards, since the dynamic flow and map can be reduced to the unit cell 
with periodic boundary conditions plus an integer value to label its position, the phase space 
can be decomposed as M. = V <g) S 1 ® Z d where V is the unit cell fundamental domain and 
d — 1 (d — 2) for a chain (bidimensional lattice) configuration. 



3.1.2 Linear stability and chaos 

The most common characterization of chaos is sensibility to initial conditions, which means 
that, given two initially infinitesimally close initial conditions, the local dynamical map will sep- 
arate them exponentially. This is characterized by the Lyapunov exponents defined as follows. 
Let Xo and Yq = Xq + 5Xq be two initially close phase space point, i.e. ||5-X"o|| ^ 1 f° r some 
norm || • ||, then 

5X t = d x &(X)\ x=Xo ■ 5X + 0{5Xl) (3.6) 
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Figure 3.2: Schematics of the collision map over a convex scatterer. A particle initially at 
X t = (q t ,v t ) collides with an obstacle boundary dD at s — qt+T(x t )- The collision occurs at 



time t + T(X ) and change the particle velocity direction according to (3.4) and (3.5). 



is the separation at time t, 5X t 
e = 5X /\\5Xq\\ is given by 



$'lo — ^'^o- The Lyapunov exponent at Xq for direction 



A(X ,e) = lim lim - In j|f *1 

<->oo pXoll-X) t \\0X \ 



(3.7) 



hence ||&Xt|| ~ e ||<5Xo||. In general dynamical systems, according to Oseledets theorem, 



(3.7) takes values from a discrete set Ai(X ) > \2(X ) > ■ ■ ■ > X r (X ) called the Lyapunov 
spectrum, whose multiplicities mi,m2, ...,m r sum up to the dimension of the phase space 



M.. Since Hamiltonian systems define symplectic systems [see (3.2)], it can be shown that the 
fundamental matrix dx^ f (X) is in fact symplectic, therefore its eigenvalues comes in pairs ±A. 
Hence, in a general Hamiltonian system with d degrees of freedom there are d independent pairs 
of Lyapunov exponents ±Aj with i = 1, . . . , d. From this follows that the sum of all Lyapunov 
exponents is null, a property satisfied by any conservative system (not necessary Hamiltonian). 
Additionally, it can be easily seen that the Lyapunov exponent parallel to the flow direction in 
phase space is null since local divergence in this direction can not be exponential; the Lyapunov 
exponent paired with the latter corresponds to the direction orthogonal to the energy surface 
H = E (this can be generalized to any constant of motion present in the system). These 
restrictions leave d— 1 not null independent Lyapunov exponents pairs. Thus, in a bidimensional 
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billiard there is only one pair ±X(X) of non-null exponents. Basically, a billiard is called chaotic 
when X(X) > for all X G M.. In case this holds in a subset of M. the dynamics is called 
mixed and otherwise is integrable. The latter occurs only for highly symmetrical billiards, the 
mixed case being the more common for an arbitrarily chosen billiard shape. 

The union of all local tangent spaces in phase space is called tangent (bundle) space and is 
where 5X lives. According to the Oseledets theorem, the Lyapunov spectrum is associated to 
a vector basis in this tangent space at each X, which can be obtained by the solution to the 
linear problem 

d x &(X) ■ ei(X) = A t (t, X) ei(&X). (3.8) 
Then, the Lyapunov exponent 

Xi(X) = lim - \n\Ai(t,X)\ (3.9) 

t— s>oo t 

is associated to the vector e.i{X). This association allows a separation of the tangent space into 
a local stable subspace, spanned by £ S (X) = {e i (X)}A l <o ! a local unstable subspace, spanned 
by £ U {X) = {ei(X)}x t>0 , and a neutral subspace £q(X) = {ei(X)}\ i=0 . Initial conditions in 
S U (X) (£ S (X)) will diverge (converge) exponentially to X. 

3.1.3 Hyperbolicity 

A stronger and more formal definition of chaotic systems is given by hyperbolicity. Given 
an invariant set A = <& l A CM, we say it is hyperbolic if for each X E A 

1. the local neutral subspace £o(X) contains only the direction of the flow, and 

2. the angle between £ S (X) and £ U (X) is always different to zero. 

A system is hyperbolic if it contains a single invariant hyperbolic subset. If the Lyapunov 
spectrum is the same for all trajectories of the invariant subset then the system is uniformly 
hyperbolic and nonuniformly hyperbolic otherwise. The latter is the most common case and 
holds in particular for hyperbolic billiards. If the flow is hyperbolic then the same holds for the 
associated map. The definition of hyperbolicity for a map is the same as before but the neutral 
direction can be ignored choosing an appropriate Poincare surface. 

In hyperbolic systems all periodic orbits have positive Lyapunov exponents which imply 
they are unstable. This is the key property to have fast correlation decay and good mixing 
properties as will be discussed bellow. Nonhyperbolic systems have stable periodic orbits which 
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live in KAM tori surrounded by unstable chaotic trajectories. The presence of these stable 
periodic orbits in phase space usually cause correlations to decay slower than in hyperbolic 
systems. 



L 




Figure 3.3: Lorentz Gas unit cell of an infinite-horizon (a) and finite-horizon (b) configuration. 
The most common form of the Sinai Billiard is given by (a). The Bunimovich Stadium (b) is 
composed of two semi-circles connected by parallel segments of length L. 



3.1.4 Sinai Billiard and Bunimovich Stadium 

Probably the most studied and well known billiard systems are the Lorentz Gas (which is 
an unfolded Sinai Billiard) and the Bunimovich Stadium. Sinai Billiards corresponds to a class 
of closed billiards composed of smooth convex scatterers (usually circles); the Lorentz Gas is 
generated by placing this unit cell periodically or randomly on the plane. The Sinai Billiard 
dynamics is completely dispersing or defocusing because an incident beam of particles gets 
defocused -the beam spreads- after colliding with the convex obstacles. It has been rigorously 
proved that Sinai Billiards are hyperbolic [Sin70j, which allows to establish strong mixing and 
ergodic properties as we discuss in the next section. 

The Lorentz Gas can be divided into two different classes depending on its geometry. When 
the obstacles are arranged in such a way that unbounded free flight trajectories can exist in 
the system, the billiard is said to have infinite-horizon. In the opposite case, it is said to have 
finite-horizon. The speed of correlations decay and hence statistical properties of hyperbolic 
billiards crucially depends on this distinction. In particular, in a finite-horizon Lorentz Gas 
the dynamics can be shown to be diffusive whereas in the infinite-horizon case the diffusion is 



anomalous [see section 3.2.2 . 
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On the other hand, the Stadium is defined as a closed billiard composed of two parallel 



line segments and two circle arcs as shown in figure 3^3 The Stadium is a mixture of focusing 
elements given by the arcs (beams of particles shrink after colliding with them) and neutral 
parallel segments where there is no dispersion. In the limit case of null length parallel segments 
the billiard reduces to a circle, which is integrable. When the circle is deformed making these 
segments sufficiently large the system becomes hyperbolic, as was proved by Bunimovich[Bun79j. 
This is explained by the fact that the focusing effect caused by the arcs is transformed into 
defocusing when the distance between them is sufficiently large. 

The focusing and defocusing mechanisms are the two basic ingredients for chaos in dynamical 
billiards. 

3.2 Statistical properties 

3.2.1 Statistical ensembles 

We now consider an ensemble of iV trajectories in phase space 

X ( t i] = e M , i = 1, . . . , N. (3.10) 

These trajectories, which form a cloud of points in phase space at each time t, do not interact 
in any way since they are different realizations of the dynamical system in question. Given an 
observable A(X), its average over the ensemble {X^} is given by 

1 N 

W'^tfE" (**')■ (311) 

i=i 

In the limit iV — > oo, the ensemble can be represented by a probability density distribution 

1 N 



i=i 



hence, the average (3.11) may be expressed as 

(A) t = [ A(X)f t (X)dX. (3.13) 

J M 
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Since probability must be conserved, from the continuity equation in phase space we obtain 

d t f t (X) = -d x ■ (F(X) f t (X)), (3.14) 

which is known as the Liouville equation. Assuming the vector field F(X) is time- independent, 
we can write its solution formally as 

f t (X) = e Lt f (X), (3.15) 

where L(f) = —dx • (F(X) /) is called the Liouville operator. In an invertible conservative 
system, the solution to ( |3.15[ ) may be expressed as 

f t (X) = f (*-*X). (3.16) 

Every phase space probability density f(X) induces a measure fi(U) = f u f(X)dX for 
U C M.. An important class of measures is induced by the stationary solutions of the Liouville 



equation (3.14), Lf(X) = 0. Invariant measures satisfy 

Ai(*-*E0 = fi(U) (3.17) 

for any phase space subset U. In general there are many invariant measures, for instance a 
density defined over any set of periodic orbits will induce an invariant measure. However, if 
the system is chaotic and we take a generic initial distribution in phase space, we might expect 
that the dynamical instabilities will produce the ensemble to be distributed over the whole 
phase space after enough time. Intuitively, the chaotic dynamic will blur any particular initial 
distribution features, and correlations between observables evaluated at the initial and long time 
densities will be null. This is a property called asymptotic stationarity, and can be formalized 
as 

\im(A) t = ! f i OD (dX)A(X), (3.18) 
Jm 

for some observable A(X) with fi^ the invariant measure. In bounded systems (with a regular 
invariant measure), the previous property is equivalent to the mixing condition, namely 

lim (A(<f> t X)B(X)) OQ = (A) 00(5)00, (3.19) 

t— >oo 

where B(X) = f (X)/f oo (X). The mixing condition implies ergodicity, which is the equiva- 
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lence of time average and phase-space average, 

lim - / A(® t X )dt = [ A(X)/i 00 (dX) (3.20) 
t Jo Jm 

for any initial condition Xq except for a set of null measure. 

In hamiltonian systems, volumes in phase space are preserved, thus the Liouville measure 
dX = d*qd*p, with (q,p) the canonical variables and / the number of degrees of freedom, is an 
invariant measure. For systems with conserved quantities such as energy, the invariant measure 
can be defined in each energy shell H = E (or equivalent) creating the microcanonical invariant 
measure. In billiards, this is given by the natural (Lebesgue) measure in Birkhoff coordinates 



X = (s, sin0) [see discussion around (3.5)] 



d 2 x = dsdsinO. (3.21) 

3.2.2 Correlation decay and diffusion 

The time-correlation function of classical observables plays an important role characterizing 
the dynamics of a system when no analytical results are known. Given two square integrable 
observables A(X) and B(X), their time-correlation function is given by 

C AB (t) = (A(&X)B(X)) - (A)(B). (3.22) 

In this section we always take averages using the relevant invariant measure, (■) = (^oo, unless 
stated differently. In systems satisfying the mixing condition, it holds that Cab(£) — > as 
t — > oo for all square integrable observables A and B. The decay rate of Cab(£) depends on how 
strongly the system mixes, i.e., it is conditional to the strength of the dynamical instability or 
chaoticity and to the smoothness of the observables. For instance, in an hyperbolic billiard with 
finite horizon (such as the Lorentz Gas) the decay of the velocity (and other smooth observables) 
is exponential. On the other hand, in systems with mixed phase space where integrable tori 
live immersed in an ergodic sea the correlation decay may be subexponential and in general 
polynomial. 

The discrete-time-correlation function can be defied replacing the flow with the collision map 



in (3.22). It is worth noting that the decay rate may be of different character for the continuos 
and discrete time correlation functions. A good example of this is the infinite horizon Lorentz 
Gas where the discrete time correlations decay exponentially as in the finite horizon case but 
the continuos time correlation decay algebraically as 1/t as a consequence of the existence of 
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unbounded free flights in the system |Ble92| . 

In systems where correlations decay sufficiently fast, the dynamical fluctuations of observ- 
ables converging to the ergodic limit are Gaussian in the sense of the central limit theorem 
(CLT). Before giving a precise definition for the CLT lets consider AA 2 = ((SU(n) — n(A)) 2 ) 
where Sa(ti) = Yl^o^-i and A n = A(<fi n Xo)- According to the discrete version of the ergodic 



theorem (3.20), in an ergodic system 



lim — SaM 

n— »oc 77, 



(A). 



(3.23) 



Then, AA 2 is the average squared of S^n) fluctuations around its ergodic limit. It is easy to 
see that, 



n-1 



AA 2 = nC AA (0) + 2 £(n - i) C AA {i). 
i=i 

Therefore, if the correlation function decays sufficiently rapid so that the surrQ 

oo 

y^CAAjn) < oo, 



(3.24) 



?1=0 



we obtain that 



where 



((S A (n) — n(A)) 2 ) = na 2 + o{n) for large n, 



a 2 = C AA (0) + 2j2C AA (n). 

n=l 



(3.25) 



(3.26) 



(3.27) 



Hence, S A (n) = n(A) + 0(y/n) for n — > oo, i.e., the average squared deviation of S A (n) from its 
mean grows as a^fn. Taking the observable as the n-th displacement A(x n ) = Ar n = r n — r n _i 
and assuming (Ar n ) = 0, we have that if the system is strongly mixing (in the sense discussed 



above) from (3.26) 



(r 2 ) = na 2 for large n, 



(3.28) 



which is the weakest characterization of deterministic diffusion in a dynamical system. We can 
identify the discrete time diffusion coefficient as D = a 2 . 



Equation (3.26), which corresponds to the Law of Large numbers in probability theory, gives 



us the average size of an observable fluctuations around its ergodic limit. A stronger property 
is given by the CLT and concerns the distribution of these fluctuations. Given an ergodic 



* According to Chernov [CheSS], the weaker condition J2^=o \C-AA(n)\ < oo also guaranties (3.261 
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dynamical system, the CLT holds if 

lim/Jx : jL= K -± <y \= \ e - z/2 dz, (3.29) 



t— >oo 



oo 



where Da is the (generalized) diffusion coefficient which can be written in continuos time as 



f 




2\ 



D A = Km - ( ( / A(&X) dt-t(A)) ) (3.30) 



o 



or, in analogy to (3.27) 



D A = / C AA (t)dt, (3.31) 

J — oo 

assuming that C AA (t) decay faster than t . When this condition does not hold the system is 



not diffusive but may be super- diffusive. Equation (3.31) is the celebrated Green- Kubo formula 
which is fundamental in non-equilibrium statistical mechanics[Dor99j. 

If we choose the observable to be the instantaneous velocity in direction x, A(X) = v x , then 



the CLT expressed in (3.29) means that the limit 



lim (3.32) 

t^oo Jdi 



exists and converge in distribution to a standard Gaussian random variable. This implies that 
[xD = Dt for large t. If the random variable 

y(s) = lim ^3 (3.33) 

is Gaussian with variance yfs for all s, then y(s) is called a Brownian process. This implies that 
particle ensembles in the billiard spread according to a diffusion equation when the billiard is 
scaled such that the mean free path and obstacles are small and the number of particles is large. 

The CLT and convergence to a Brownian motion has been proved for general finite-horizon 
hyperbolic billiards by Bunimovich, Sinai and Chernov[BSC91J. An analogous result was 
obtained for the infinite-horizon Lorentz Gas first by Bleher [Ble92j and then formalized by 
Szasz [5V07j : in this case a non- normal CLT was proved with 

t^OO y/DtlOgt 
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the correct stationary random process. This implies diffusion is anomalous, with (x 2 ) ~ t\ogt. 

3.3 Machta-Zwanzig approximation 

In strongly chaotic periodic billiards with a space configuration such that the escape time 
from a unit cell is longer than its ergodic time, i.e. the time mixing takes place within a unit 
cell, the motion of particles in the periodic lattice is well approximated as a (symmetric) random 
walk. Indeed, if a typical trajectory is trapped within a unit cell for a time sufficiently long so 
that mixing takes place, then the decay of correlations imply the particle will lose memory of its 
initial state; therefore, its exit direction (into a neighboring unit cell) will be effectively random. 
We note that this excludes unit cells with infinite-horizon geometries. Expressed formally, given 
an initial phase space point X = (q , ^o? Io) £ ® S 1 <8> when the mentioned conditions 
hold the unit cell index (winding number) dynamics I n follows a Markov process, i.e., 

P r (-^n+l | Im In-li ■ ■ ■ : Iqi ■ • ■) — P r (-^n+l I In)- (3.35) 

Under this approximation is possible to give an analytical expression for the diffusion coefficient 
as we discuss below. 

The idea to obtain the billiard diffusion coefficient in this way was first considered by Machta 
and Zwanzig for the particular case of the Lorentz Gas in the high density regime|MZ83j. In 
this seccion we review their argument for a general one-dimensional billiard chain assumed to 
be strongly mixing and to possess narrow exits compared to the dimensions of the cavity. In an 
ergodic planar billiard, the average residence time in a unit cell T>q is given by |Che97j 



\dv a 



(3.36) 



where |X> | is the fundamental domain area and \&Dq\ is the exit openings total length [for 



instance, consider the unit cell given by figure 3.3 (b) where the dashed lines correspond to dV ]. 



On the other hand, it is well known that for a random walk in an isotropic]^] unidimensional lattice 
with period a the diffusion coefficient is D = a 2 jr. Therefore, we obtain the Machta-Zwanzig 
diffusion coefficient, 

D MZ = (3.37) 

7T \V \ 



^This means that the probability of a particle to exit through the right and left leads are the same and equal 
to 1/2 independent of its prior direction. 
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with Z the unit cell length. We note that (3.36), usually called Santalo's formula, is an exact 



result so the only assumption to obtain (3.37) was the random walk approximation. As we have 



argued, this condition holds when the particles are trapped for a long time in a unit cell beacuse 
this means they collide many times with the obstacles; hence, since the billiard is assumed 
strongly chaotic this implies a fast correlation decay and effective markovian dynamics for the 
winding number J n leading to the random walk approximation. Consequently, assuming strong 



chaotic dynamics in the unit cell, the validity of (3.37) rests only on the mean free path length r, 



since the average number of collisions of a trajectory before exiting a cavity is given by t/t c . An 



analytical expression for r c is given by (3.36) replacing |<9X>o| with |0Pq|> the obstacles boundary 
length in the unit cell. Then, we have that the Machta-Zwanzig approximation for the diffusion 



coefficient (3.37) holds in the limit r c <C r, that is, when 



\dv \ « \dv \. 



(3.38) 



3.4 Cosine billiard 

In this section we define the periodic cosine billiard which will be used in the next chapters 
to numerically test some of our classical and quantum (wave) results; therefore we sometimes 
also refer to it as the cosine waveguide. The periodic cosine billiard is composed of a one- 
dimensional chain of hard-wall cavities with cosine shaped boundaries. We employ this system 
because our numerical method for the quantum scattering problem -which is very efficient in the 
semi-classical limit- requires a smooth waveguide with connected boundaries, thus the Lorentz 
Gas or Bunimovich Stadium chain cannot be employed. Even though the quantum numerical 
method is highly efficient for our purposes, the cosine billiard has the disadvantage to be much 
harder (numerically expensive) to solve classically that the aforementioned systems. In addition, 
there are no analytical proofs of hyperbolicity and ergodicity for the cosine billiard, so we will 
have to investigate them numerically. In spite of these contretemps, the cosine billiard has 
been used before in the quantum chaos literaturejIMRK961 IHKLOOl IMBLAicvP02l IMSVL08| 
because of the possibility of easily changing its dynamic from a mixed phase space to (apparent) 
full chaos with the tuning of the cosine amplitude. 

Let (x, y) be the coordinates in configuration space. We recall we have set the particle's 
speed to v — 1 unless otherwise stated; for general velocity v the diffusion coefficient simple 
scales asD-> Dv. We define the unit cell as the region enclosed by h\{x) < y < h,2(x) for each 
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Figure 3.4: Cosine billiard chain with five unit cells connected to two plane leads. The unit cell 



boundaries are defined in (3.39)-(3.40) as a function of the amplitudes Ai and A 2 shown in the 
figure. 



x G [—1, 1], where 



A\ 

h\(x) = — [1 + cos (rex)] and 
A 2 

h 2 (x) = A\ + — [1 + cos (nx)} , 



(3.39) 
(3.40) 



Hence, h\(x) and h 2 (x) define the lower and upper waveguide boundaries, respectively. With 
the unit cell defined in this way our cosine billiard always has finite horizon, i.e., it does not 
allow unbounded collision-free trajectories for any values of A\ > and A 2 > 0. It is trivial 
to see this is guaranteed by the fact that hi(0) = h 2 (±l) = Ax. As explained in section 3.1.4 



besides strong mixing, the finite-horizon property is fundamental to obtain normal diffusive 
dynamics. 

Note that the unit cell mirror symmetry x — > — x is not relevant for the classical transport 
properties of the billiard but makes the numerical solution of the quantum scattering problem 
faster (the transmission and reflexion matrices t and r are the same in both direction). However, 
this induces an anti-unitary symmetry in the quantum Hamiltonian which plays a role in the 
statistical and transport properties of the waveguide as we discuss in chapter [5} 

The character of the cosine billiard classical dynamic can be mixed or predominantly chaotic, 
depending on the parameters A\ and A 2 . In order to assess the dynamical properties of this 
system we use basically two tools, namely the unit cell Poincare sections 



{(y, v y ) : WneN, \x\ = 1 + nZ} 



S y = {(s,v s ) : Vx(s) e [-1,1], y = h 2 [x{s)}} 



(3.41) 
(3.42) 
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where Z = 2 is the waveguide period and (s, v s ) denote the Birkhoff coordinates, and the 
instantaneous velocity autocorrelation function 

C v (t) = (v x (t)v x (0)). (3.43) 

We note that S x is the stroboscopic Poincare section on the unit cell connecting border and 
S y is the surface of section on the unit cell upper wall parametrized with Birkhoff coordinates. 
Transport properties of the cosine billiard depend mostly on its dynamic over S x , because when 
this map is sufficiently chaotic almost all stable quasi-periodic trajectories on S y are trapped 
within the unit cell and therefore do not contribute to transport. 

We will consider two ranges of parameters. The first is given by Ax = 1 and A 2 G [0.3, 0.6] 
where the system is predominantly chaotic but a few small tori are observed in S y for some values 
of A 2 [see figure [3~5| . Still, in all cases tested in this range, the connected ergodic component 
in S x makes up 95% or more of the section area and C v (t) decays exponentially [see figure 



3.6 . This implies normal diffusive dynamics as can be seen by looking at the scaled asymptotic 



displacement 

x = hm — , (3.44) 

t-KjC y/t 

which is distributed as a Gaussian random variable with zero mean and variance D = {x 2 ) 



[see figure 3.7 . Since the velocity of autocorrelation decay is exponential, the convergence of x 
distribution to its asymptotic Gaussian limit is quite fast. An alternative method to support 
our conclusion that x/\/Di converges to a standard Normal distribution and that at least its 
first two moments also converge, is given by calculating Mi = (\xt\). It is easy to see that if 
this is the case then 

^-M 2 = {x 2 } = D (3.45) 



2 



see figure 3.8 . Note that, for instance, an anomalous diffusive system with ergodic phase-space 



but infinite horizon trajectories may exhibit convergence in distribution to a Gaussian for an 



appropriate normalized displacement but would fail to satisfy relation (3.45) between its first 
and second moments for any finite time approximation of x [AHO02J. 

As we have seen for the parameter range A\ — 1 and A 2 € [0.3, 0.6], although there are small 
tori observable in the section S y , the billiard dynamic is strongly chaotic in S x and effectively 
diffusive as we have seen from the fast decay of correlations C v (t) and strong convergence (at 
least up to the second moment) of the scaled displacement x. In chapter [5] we will use this 
configuration to test a semi-classical result regarding the average number of propagating Bloch 
modes in a diffusive waveguide and show excelent agreement with the analytical result. However, 
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this range of parameters has the disadvantage of not fulfilling condition (3.38) for the Machta- 
Zwanzig approximation. In fact, it is easy to see that \dT>o\ = 3A\ + A2 and \&Dq\ = A\, thus 
in the best case \dT> \ = 3.6\dV \. 

The second parameters range we consider is given by A\ = 0.5 and A2 G [2.5, 4.5], which 
is closer to the Machta-Zwanzig limit since in this case |(9Po|/|<9^o| moves between 8 and 12. 
With this configuration the system satisfies all the properties previously discussed, i.e., it shows 



fast decay of correlations and exhibits a diffusive dynamic [see figures 3.5 and 3.6 
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Figure 3.5: Density plots of the Poincare sections S x (right column) and S y (left column) of the 
cosine billiard. Figures (a) correspond to a configuration with parameters A\ = 0.5, A 2 = 3.0 
where we observed strong chaotic dynamics; in this case there are no visible tori. In (b) we 
show a configuration with A\ = 1.0, A 2 = 0.45, where the dynamics is also strongly diffusive 
but a small tori is present in S x related to families of integrable stable orbits trapped in the 
unit cell; this is not relevant for transport since the stroboscopic map is strongly chaotic as the 
density plot in S y reveals. 
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Figure 3.6: Velocity autocorrelation funcion C v (t) = (v x (t)v x (0)) for an ensemble of 10 5 particles. 
In figure (a) we plot this function for the three configurations with Ai = 1.0, A 2 = 0.6, 0.45, 0.3 
(full, dashed and dotted lines). In figure (b) we show the same function for Ai = 0.5, A 2 = 
2.5, 3.5, 4.5 (full, dashed and dotted lines). The decay to noise level is quite fast in all cases; in 
(a) at around t ~ 20 (which is approximately 25 collisions) correlations have decayed to noise 
level whereas in (b) the same happens for t ~ 30 (which is approximately 9 collisions). 
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Figure 3.7: Histogram of the normalized displacement x t for a cosine waveguide with A\ = 1.0, 
A 2 = 0.3 at time t = 50000 using an ensemble of 10 5 initial conditions. The best Gaussian fit 
is plotted over it. The inset show the same histogram (dots) in log-scale where convergence to 
a Gaussian (full line) is more evident to be achieved even deep in the tails. For the A\ = 0.5 
configurations of the cosine billiard we are considering a similarly good converge of x t is observed. 
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Figure 3.8: Finite time approximations of the diffusion coefficient D(t) = {xf)/t for the A 1 = 
1.0, A2 = 0.3 cosine waveguide calculated from the variance (second moment) of the sample's 
displacements (black line) and from the first moment using §(|£t|) 2 = (xf), as expected for 
Gaussians, (dashed line) for an ensemble of fO 5 initial conditions. To obtain the diffusion 
coefficient D we average D(t) for times after the transient relaxation regime. The diffusion 
coefficients obtained in this way from the first and second moments of x t are the same up to 
standard error. 
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Chapter 4 

Random matrix theory 



The theory of random matrices was first formally developed in the nineteen-sixties mainly 
by Wigner, Dyson, Mehta and Gaudin, motivated by understanding the energy spectra statis- 
tics of heavy nuclei. Later in the same decade, these techniques were applied to small metal 
particles to study their microwave absorption properties. A brief account of the history of these 
developments as well as a complete treatment of the theory's current status can be found in the 
book by Metha |Meh04| . 

Random matrix theory (RMT) experienced a revival of interest in recent decades and has 
found many applications in disciplines outside physics, for instance in statistics, finance, en- 
gineering and number theory. In physics, its domain of use shifted from nuclear physics to 
quantum chaotic systems since it was noted by Bohigas, Giannoni and Schmit|BGS84j that the 
Wigner-Dyson ensembles of hermitian matrices applied generically to describe statistical prop- 
erties of closed chaotic systems (see also |Haa01] and |Gut90j ). This discovery was followed by 
the work of Altshuler and Shklovskii|AS86j on the universal conductance fluctuations of disor- 
dered systems, which led to the development of a random matrix theory of quantum transport. 
A nice review on this subject is given by Beenakker in |Bee97j . 

In this chapter, we present some known properties and results of RMT which are relevant 



for our work. First, in section 4.1 we discuss the connection between a quantum chaotic 
systems and random matrices. Then, in the following sections, we present the Wigner-Dyson 
and Circular ensembles which describe the Hamiltonians and Scattering matrices of chaotic 



systems, respectively. Finally, in section 4.3.4 we define the RMT periodic chain model which 
we will use in the next chapters. 



4.1 Connection to chaotic systems 

Random matrix theory deals with the statistical properties of large matrices with randomly 
distributed elements. The basic assumptions for RMT to be a good description of a given 
physical system are that its evolution dynamics maintains phase coherence, it is chaotic and 
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it is big so that the number of energy levels (or degrees of freedom) is sufficiently large for 
a statistical description to make sense. Hence, for low-dimensional chaotic systems such as 
billiards, the connection between their quantum statistical properties and RMT is achieved in 
the semi-classical limit. 

Let us consider a quantum system whose energy spectrum is E\ < Ei < ■ ■ ■ < En for 
a fixed integer N which is large but finite. One of the first spectral statistics to be studied 
semiclassically was the energy level nearest-neighbor spacing P(S), which is the pdf of the 
unfolded [HaaOl] spectrum nearest-neighbor spacings, namely 

Si = Ei ~l i+ \ (4-1) 
AE K J 

where AE = 1/pE is the mean level spacing with pE the mean local density of states [see 
equation ( 5.9[ )]. Generic classically integrable systems with more than two degrees of freedom 



have a Poissonian level spacing distribution P(S) = e~ s and display no correlations, hence 
energy levels tend to cluster together without any repulsion |BT77j . This implies that level 
crossings are not avoided when a parameter in the Hamiltonian is changed. On the other hand, 
classically chaotic (non-integrable) systems display different levels of repulsion depending on 
their symmetries; there are three universality classes with repulsion given by 

P(S) ~ for S -+ 0, (4.2) 

with = 1,2 or 4. For systems possessing an anti-unitary symmetry (such as time reversal) 
= 1, for systems without such symmetries = 2 and for anti-unitary symmetric systems with 
broken spin-rotation symmetry /3 — 4. It was noted |BGS84] that the full distribution P(S) ob- 
served in chaotic systems corresponded for all S to the nearest-neighbor spacings observed in the 
eigenvalues of large random hermitian matrices taken from the appropriate gaussian ensembles, 
called the Wigner-Dyson ensembles. Thus, it was conjectured by Bohigas, Giannoni and Schmit 
(BGS) that a system with chaotic phase-space dynamics is expected to show universal quantum 
level statistics in the semiclassical limit, consistent with the predictions of these ensembles, with 
the only remaining relevant physical parameter being the mean level spacing AE. 

The main tool to understand this connection between Wigner-Dyson RMT and chaotic 
systems has been the Gutzwiller trace formula jGut90], which was recently used in the first 
semiclassical proof of the BGS conjecture, given by Hakee et al. in |MHB+04l IHMA+07] . 
Gutzwiller trace formula gives a semiclassical expression for the density of energy states as a 
sum over the periodic orbits of a chaotic system, p(E) = A~ t e lS " , l r \ where S-y and are the 
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action and stability amplitude of the periodic orbit 7. After unfolding the spectra of a chaotic 
system, all statistical properties of p(E) are reproduced by the analogous function in an ensemble 
of random Hamiltonians with the appropriate symmetries. In disordered systems, the physical 
ensemble represented by RMT is given by the set of allowed disorder realizations. In contrast, in 
chaotic systems where disorder does not play a role the physical ensemble being caricatured by 
RMT can be taken as a uniform set of energy realizations in the interval [E — dE/2, E + dE/2], 
with AE <C dE <C E such that p(E) is approximately constant in this interval. This means 
that dE is classically small but large enough to contains many quantum levels. We call this the 
semiclassical ensemble. 

On the other hand, the RMT of quantum transport is based on Dyson circular ensembles of 
scattering matrices, which can be defined by a principle of maximum entropy on a set of unitary 
matrices with the appropriate symmetries imposed (unitary, orthogonal or symplectic). First, 
Bliimel and Smilansky[ BS9Q ] showed semiclassically that the distribution of scattering phase 
shifts of a finite-horizon chaotic cavity was well described by the circular ensembles and then, 
Jalabert, Pichard and Beenakker[JPB94] showed that this RMT agreed with the semiclassical 
conductance in a quantum dot calculated by Baranger, Jalabert and Stone[BJS93j. This theory 
explains two of the most salient features of quantum transport of chaotic cavities, namely 
the weak localization effect and universal conductance fluctuations. The transport properties 
in a disordered wire geometry was first studied by Dor okhov |Dor 82] and Mello, Pereyra and 
Kumar |PAPN88] . They developed a differential equation for the evolution of the transmission 
eigenvalues pdf as a function of the wire length known as the Dorokhov-Mello-Pereyra-Kumar 
(DMPK) equation, which assumed a wire composed of small weakly scattering segments taken 
from the circular ensembles. This equation accounted for the ohmic and localized regimes 
observed in a disordered wire. Another method used in quantum transport is the non-linear a 
model, which allows us to compute the n-point correlation functions of the wire eigenvalues; this 
method has been proved equivalent [BF96J to the DMPK equation but provides less information 
since only give access to correlation functions. A different approach was later pursued by 
Iida, Weidenmuller and Zuk |IWZ90at IIWZ 90bj whom computed the conductance average and 
variance of a (discrete) chain of random quantum dots taken from the circular ensembles. They 
studied the relevance of the coupling of the quantum dots chain to the external leads and 
showed that, assuming orthogonal Wigner-Dyson statistics for the quantum dots hamiltonians, 
was sufficient to obtain universal conductance fluctuations. 

Later, Simons and Altshuler |SA93al ISA93bj considered the statistical properties of the spec- 
trum Ei(X) when a chaotic system is subject to an external perturbation X, for instance in 
the form of a Aharonov-Bohm magnetic flux or a background potential. In this case, it is 
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found that the statistical properties of Ei(X) are universal after appropriate scaling depen- 
dent on two physical parameters: the mean level spacing AE and the generalized conductance 
C(0) = (dxEf)/(AE) 2 . This allows us to examine the statistical properties of the energy bands 
in periodic systems[M CAJ94] . where X is the Bloch phase which acts as an Aharonov-Bohm 
flux over the unit cell defined on a torus topology (periodic boundary conditions). In section 



5.1.3 we will employ Simons and Altshuler universal correlation function to obtain the average 



number of propagating Bloch modes in a periodic chaotic system. 

4.2 Wigner-Dyson ensembles 

4.2.1 Definitions 

There are three different Wigner-Dyson ensembles of hermitian matrices, namely the Gaus- 
sian Orthogonal Ensemble (GOE), Gaussian Unitary Ensemble (GUE) and Gaussian Symplectic 
Ensemble (GSE). They receive their names because they contain random gaussian elements and 
are invariant with respect to orthogonal, unitary and symplectic transformation, respectively. A 
detailed treatment of this subject can be found in |Meh04] : in the present section we summarize 
the relevant results. 

A generic physical system without time-reversal invariance possess a Hamiltonian not re- 
stricted by any symmetry but being hermitian. Given an hermitian matrix H = H . the 
transformed matrix 

H' = U X HU (4.3) 

is also hermitian for any unitary matrix U G U(N). Hence, the GUE ensemble is defined over 
the space of hermitian matrices with a pdf P(H) satisfying 

P(H') dW = P(H) dH (4.4) 



for any unitary matrix U with H' given by (4.3) and where the volume element is given by 



dH = HdH$l[dH§ (4.5) 

k<j k<j 

with Hjfj and , the real and imaginary parts of Hkj = (H)kj, statistically independent 
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variables. It follows from this definition that most general form of the GUE pdf is 

P(H) = cexp (-a Ti H 2 + b Tr H) (4.6) 

with a > 0, b and c real. The constant b is usually taken as zero since it accounts to a simple 
change of the spectrum Ej mean value. 

On the other hand, if a (spinless) system possess time-reversal symmetry, its eigenfunctions 
are real, hence there exists a non-diagonal basis where the Hamiltonian is a real matrix. Thus, 
in this case the Hamiltonian can be taken as a real symmetric matrix. The GOE ensemble 
is defined over the set of these matrices requiring equivalent invariance properties than for the 



GUE case, namely that (4.4) holds with H' = W l HW for any orthogonal matrix W G O(N). 



The GOE volume element is similar to (4.5) without the imaginary parts product. From this 



conditions follows that the GOE pdf have the same form as (4.6). 

The third universality class corresponds to systems with time-reversal symmetry but without 
spin-rotation symmetry. It can be shown that in this case the system Hamiltonian can be chosen 
as a hermitian self-duaQmatrix H = H R . The set of matrices with these symmetries is invariant 
to symplecticjT] transformations, hence the GSE is defined imposing this requirement in a similar 
fashion to what we have done for the GUE and CUE. 

Originally, Wigner and Dyson studied an ensemble of hermitian matrices H with P{H) = 
c exp(— a Tr V(H)) with V(-) a potential function. In case V(H) ~ H 2 we recover the Gaussian 
ensembles; however, in the limit of large matrices N — > oo the spectral correlations become 
independent of V(-) sufficiently far away from the spectrum edges. This is the heart of the 
universality of spectral correlations underlying BGS conjecture. 



4.2.2 Eigenvalue and levels spacing distributions 

Once we have obtained the pdf for the Gaussian ensembles, the next question is to determine 
the eigenvalues {E n } pdf. This is accomplished by finding the Jacobian Jg that allows us to 
write the volume element dH as a function of E n and the eigenvectors matrix. Thus, we have 
that 

TV 

dH = Jpdfip(U)Y[dEi, (4.7) 

i=i 

*This means that the Hamiltonian is a quaternion matrix and there exists a quaternion basis where it is 
quaternion real. 

^Here we refer to the symplectic group Sp(iV) which is the group of invertible quaternionic matrices preserving 
the standard hermitian norm. 
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where /3 = 1,2,4 for the unitary, orthogonal and symplectic ensembles, U is in the relevant sym- 
metry group and dfipiU) is the Haar measur^jof this group. The Jacobian is given by |HaaOT] 

N 

M{En}) = l[\E t -E J f, (4.8) 



with (3 = 1,2,4 defined previously. Since Tr H 2 = ^2 i=1 Ef, P{H) given by (4.6) does not 



depend on the eigenvectors of H. Therefore, the eigenvalue distribution takes the form 

P{{E n }) =cY[\Ei- Etf exp l-a^tf > ( 4 - 9 ) 

i<j V t=l / 

which can be recasted in the form of a Boltzmann-Gibbs distribution as 

N N 



P({E n }) = c exp 



-P [ J^uiE^E^+aJ^Ei 

. i<3 i=l 



(4.10) 



with a two-point potential energy u(E,E') = — ki\E — E'\ and (5 playing the role of inverse 
temperature. Hence, the energy eigenvalue dynamics when a system parameter is changed 
is expected to be analogous to a unidimensional Coulomb gas since the logarithmic repulsion 
u(E,E') is equivalent to the Coulomb interaction of two line charges. Although the eigenvalue 
dynamics turns out to be integrable, their statistical equilibrium properties are reproduced by 
RMT. 

The marginal probability density function 



p(E) = JdE 2 ... dE N P(E, E 2 , . . . , E N ) (4.11) 



is also of particular interest. Setting the normalization (HijHji) = 1/4JV (fixing the constant a 
in ( 4.6| )) in order to avoid divergences in the limit iV — > oo, one obtains the well-known Wigner 



semicircle law 

, s f (2/tt) Vl-E 2 for \E\ < 1 , 

m = \ y } . ' ' (4-12) 

[ for \E\ > 1 

valid for the unitary, orthogonal and symplectic ensembles. 

Finally, using P({E n }) we can compute the nearest-neighbor spacing pdf. For the most 

■^Given a topological group G, the left (right) Haar measure is the unique measure over subsets of G that is 
invariant to left (right) multiplication. For compact groups (such as U(iV), O(N) and Sp(iV)) the right and left 
Haar measures are the same. 
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trackable case of dimension N = 2 is possible to obtain a simple closed form solution, namely 



717 



rS/2)e-^/ 4 orthogonal 
P(S) = { (32S 2 /n 2 )e- is2 / 7r unitary (4.13) 

(2 18 S73 6 7r 3 )e- 6452 / 9 " symplectic. 



The distribution P(S) in the limit N — > oo is quite close to (4.13) so for most practical purposes 
this is an adequate approximation. In particular, it gives the correct level repulsion (4.2) for 
the three universality classes. 



4.3 Circular ensembles 

4.3.1 Definitions 

The circular ensembles are probability spaces defined over unitary matrices satisfying a 
maximum entropy principle, i.e. their probability measures are uniform in the appropriate 
parametrization sets. In the same fashion as with the Wigner-Dyson Gaussian ensembles there 
are three classes of circular ensembles, namely the Circular Unitary Ensemble (CUE), Circu- 
lar Orthogonal Ensemble (COE) and Circular Symplectic Ensemble (CSE). These ensembles 
were originally developed by Dyson |Dys62| as a more manageable alternative to the Gaussian 
ensembles but later a microscopic justification was given by Brouwer |Bro95] relating them to 
the Gaussian ensembles; given a quantum dot whose Hamiltonian is taken from GUE, GOE or 
GSE, its scattering matrix belongs to CUE, COE or CSE, respectively. 

The scattering matrix of a system without time-reversal symmetry is restricted only to be 
unitary. Since any unitary matrix S e U(iV) can be written as S = UV, with U and V also 
unitary, we define a neighborhood dS in the set of unitary matrices as 

S + dS = U(l + idM)V, (4.14) 

where dM is an infinitesimal hermitian matrix with elements dM^ = dM^y + idM>y and we 
note that here N is the total number of scattering channels in the system (for instance, in a 
two-leads waveguide with M\ propagating channels on lead 1 and M2 on lead 2, N = Mi + M2). 
The real and imaginary parts of dM are allowed to vary independently over small intervals 
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dfifj , d/j,fj , so the volume element measure is 

N N 

yUU (rf5)=n^g ) ri rf 4 ) ' ( 4 - is ) 

i<j i<j 

independent of the choice of U and V. Since we require a uniform probability density function 
Pv(S) over the unitary group U(N), 

Pv(S)dS = -J— Mu (dfif) (4.16) 

with the normalization constant Vu(iv) given by the measure of the complete group. Dyson 



defined the CUE ensemble by property (4.14), namely requiring an ensemble pdf invariant to 
every automorphism S — > U'SV with U',V G XJ(N), which implies that fi\j is the Haar 
measure of the unitary group. 

On the other hand, for a system possessing time-reversal invariance and spin-rotation sym- 
metry, the scattering matrix must be a unitary symmetric matrix S, which can always be 
written as 

S = U T U (4.17) 
with U G U(N). Thus, in this case the infinitesimal neighborhood is defined by 

S + dS = U T (1 + idM)U, (4.18) 

where dM is a real symmetric matrix with elements dMy, % < j, varying independently in an 
small interval dfaj. Then, the measure of this differential is 

N 

fio{dS) = Y[dfHj. (4.19) 

Hence, we define the COE probability density Po(S) to be uniform over the set of unitary 
symmetric matrices, this is 

Po(S)dS = -J—fMQ^dS), (4.20) 

t / Us(AT) 

where Vu s (at) is the normalization constant. Dyson proved that the COE ensemble is uniquely 



defined by property (4.18), this is by demanding an ensemble pdf invariant to every automor- 
phism S — > V T SV with V G XJ(N). Note that the set of unitary symmetric matrices do not 
form a group, thus, in particular do not have an associated Haar measure. However, using 
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decomposition (4.17) it is possible to obtain fio as an induced measure from n\j. 

Finally, for time-reversal invariant systems but with broken spin-rotation symmetry, the 
scattering matrix is a unitary self-dual quaternion matrix. Any such matrix can be decomposed 
as 

S = U R U (4.21) 

with U e U(iV) and U R its dual. In this case Dyson proved that the CSE ensemble is defined 
uniquely by the property of being invariant to every automorphism S — > W R SW where W is 
a unitary quaternion matrix. The CSE measure can be obtained in a similar way to fio, namely 



as the induced measure from nu and decomposition (4.21). 



4.3.2 Transmission eigenvalues distribution 

Since the circular ensembles pdfs are uniform, it is clear that when transforming to its 
diagonal form representation the only non-trivial term is the Jacobian J related to the eigenvalue 
variables. Given a scattering matrix S from one of the circular ensembles, we can write dS = 
J HE{dU) Y\i d<pi, with (f>i e [0, 27r) the phase-shift (i.e. the eigenvalue phases), U its eigenvectors 
matrix (which is orthogonal, unitary or symplectic, depending on the ensemble considered) and 
\1e the respective unique ensemble measure. The Jacobian, 

N 

J = JJ|e^ ! -e^\ p , (4.22) 
was also computed by Dyson and is basically the joint eigenvalue pdf, P({0/}) = c J, with c a 



normalization constant. The repulsion exponent f3 in (4.22) is the same as for the Gaussian case, 
namely /3 = 1,2,4 for the Circular Unitary, Orthogonal and Symplectic Ensembles, respectively. 
As we have seen in chapter [2j the S matrix eigenstates mix incident and outgoing waves [see 



equation (2.8)] and its eigenvalues are not directly related to transport properties. This is the 



reason why we have introduced the polar decomposition in order to relate the scattering matrix 



to the transmission eigenvalues T« [see section (2.1.3)]. Therefore, from a transport perspective 



it is more natural to consider the ensemble density as a function of these variables; in other 
words, the Jacobian associated to this representation must be computed. The corresponding 
measure in this representation was obtained by Baranger and Mello[BM94j and is given by 



N N 



i<j 
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with U n the four unitary matrices in (2.28) which in addition are symmetric or self-dual for 
COE or CSE, respectively. 



4.3.3 Quantum dot and DMPK equation 



As we have mentioned in section |4.1[ the Circular Ensembles where first employed in the 
context of quantum transport in order to model disorder cavities (quantum dots) and disordered 
wires. In the case of quantum dots, the scattering matrix is directly taken from the appropriate 
ensemble depending on the desired symmetries of the system. Then, all transport information 



is encoded in P(S) given by (4.23), in particular the conductance g = Yli^i- Since the con- 
ductance can also be written as g = ^2 nm \t n m\ 2 , an important quantity is given by (|t„ m | 2 ) for 
the Circular ensembles. In fact, from the construction of the Circular ensembles is easy to show 
that the squared average of the scattering matrix elements^ is given by 



( I S nm | ) 



1 

2N 

1 + 2N 



for CUE 



for COE. 



(4.24) 



In the CUE case, transmission and reflexion probabilities are all equal, independent of the 
incoming and outgoing channel; this is the quantum analog to the Machta-Zwanzig assumption 



of an effective random walk dynamics [see section 3.3 . On the other hand, (4.24) shows explicitly 
that in the COE ensemble scattering from mode n back to mode n is twice as probable as from 
mode n to a different channel, which is a consequence of time-reversal symmetry. 



From (|4.24|), we obtain the average conductance in a quantum dot with two identical contacts 

(g) 



with N channels, 



N 
~2 



1 + 



1 



4 



1 



This means that 



(g) 



COE 



(g) 



CUE 



1 

4' 



(4.25) 



(4.26) 



implying that the conductance in time-reversal symmetric systems is lower than when this 
symmetry is broken. This effect, called weak localization, is explained semiclassically by the 
enhancement of reflection probability due to constructive interference of time-reversed trajecto- 
ries [Dat95j. Weak localization is a purely quantum effect; classical transmission and reflection 



§Here we assume a quantum dot connected to two identical leads with N propagating modes. The general 
case with differents leads is also easy to obtain[ Bee97j . 
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probabilities off a quantum dot are equal [see section 3.3| . 

Another well known universal feature observed in quantum dots are the universal conduc- 
tance fluctuations: for two equal contacts with N channels the conductance variance is 

Va%] = ^/T\ (4.27) 

independent of N and hence of (g). 

On the other hand, the standard model to describe a disordered wire is the DMPK equation. 
This is a partial differential equation for the transmission eigenvalues pdf as a function of the 
wire length. Its derivation starts from a finite wire segment with transmission spectrum T n 
which is composed with a new segment of infinitesimal length. Then, the compounded wire 



transmission matrix is obtained using [see equation (2.13) 



t 2 = t 1 (l-r r 1 )- 1 t , (4.28) 

where the index i — 0, 1,2 denote the S matrix components of the original, infinitesimal and 
resultant wires, respectively. It is assumed that the conductor is weakly disordered (with an 
order parameter given by l/L <C 1, where L Q and I are the infinitesimal segment length and 
mean free path) so that we can write the resulting transmission eigenvalue perturbatively as 
T n + 5T n . In order to obtain ST n , a second order perturbative expansion of t 2 t\ is performed. 



Then, the transmission and reflection matrices are written in polar form [see equation (2.28)] 
with the matrices U, V, U', V taken from the desired circular ensemble. This allows us to 
compute the first two moments of ST n , which imply a Fokker-Planck equation for T n with the 
wire length playing the role of time parameter. For a detailed derivation of the DMPK equation 
see the book by Mello and Kumar |MK04| . 

The DMPK equation shows that in a disordered wire there are three regimes as the systems 
length L is increased, namely a ballistic regime for L/l <C 1, then a metallic regime for 1 <C 
L/l <C iV and finally a localized regime for N <C L/l. In the metallic regime, the average 
conductance decays linearly with L displaying ohmic scaling (g) = N/(s + 1) with s = L/l. In 
addition, in this regime the conductance exhibits a phenomenology similar to the quantum dot, 
showing weak localization and universal conductance fluctuations with 

5g = -\ (4.29) 
Va%] = ^/T\ (4.30) 
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where 5g = (g)coE — (g)cvE- Note, however, that this result differs slightly from the quantum 



dot [equations (4.26) and (4.27)], which is a consequence of correlations in the DMPK equation 



solution being not geometric as for the quantum dot in (4.23). 



The localized regime is characterized by the exponential decay of the conductance (g) = 
exp (— L/2£) with £ ~ N/31, which is also known as Anderson localization |AALR79] . This 
effect is again a consequence of wave interference but in this case not restricted to time-reversed 
trajectory pairs; multi-scattering in a long wire can produce strong inference among different 
ray paths causing the absence of transmission. In this regime, the conductance pdf is log-normal 
with Var[log£f] = -2 (logs) = 4L/£. 



4.3.4 Periodic chain of quantum dots 

The periodic chain model we use in the next chapters is constructed by joining together a 
series of identical quantum dots defined by a scattering matrix taken from the CUE or COE 
ensembles. The scattering matrix for arbitrary chain length can be obtained using composition 



rules (2.13), thus each quantum dot (circular) ensemble induces a periodic chain ensemble. Note 
that this ensemble of scattering matrices is different from the Circular ensembles because, as we 



have seen in section 2.2.2 in the long chain limit transmission eigenvalues Tj associated to non- 
propagating Bloch modes converge exponentially to zero whereas Tj linked to propagating Bloch 



modes are 0(1). This implies a distribution P(S) different to (4.23) and, in fact, dependent on 



the chain length. In chapter |6l we will study the conductance pdf of this model as a function of 



its length [see figure 6.11 where some of the particular features of the periodic chain ensemble 
become evident. 
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Chapter 5 

Number of propagating Bloch modes 



As we have shown in 2.2.1 the number of propagating Bloch modes Nb plays an impor- 
tant role for the transport properties of a long periodic chain. The asymptotic semiclassi- 
cal behavior of (Nb) depends on the phase-space dynamics of the unit cell, growing linearly 
with the wavenumber k in systems with a non-null measure of ballistic trajectories and go- 
ing like ~ V Dk in diffusive systems (with D the associated classical diffusion coefficient) as 
was shown by Faure [FauQ2j . We review this result in section 5.1, where we also present nu- 
merical evidence of its validity in a physically realistic waveguide. In addition, we show an 
alternative derivation of (Nb) using Simons and Altshuler's universal parametric autocorrela- 
tion function |SA93at ISA93bj . Finally, in section 5.2 we obtain an exact expression for (Nb) 
in a RMT periodic chain using Faure's result and the Machta-Zwanzig approximation for the 
diffusion coefficient. 



5.1 Nb in the semiclassical limit h — > oo 

5.1.1 Faure's result 

We start reviewing the main result obtained by Faure in |Fau02] for (Nb) in general Hamilto- 
nian quasi-one-dimensional periodic systems and then obtain the particular result for a billiard. 
Let us consider the stationary Schrodinger equation in coordinates (x,y), 

h 2 V 2 \ 

+ V(x,y))if)(x,y) = Eip(x,y), (5.1) 



2m 

with V(x,y) the potential assumed to be Z-periodic in x, i.e. V(x + Z,y) = V(x,y), and 
constrictive in y, this is, 

V(x, y -> ±oo) = oo. (5.2) 



As we have recalled in section 2.1.2, since the Hamiltonian is periodic we known from Bloch 



theorem that the solutions of (5.1) are given by E n (0), ip n ,e(x,y) with n G N and 6 £ [0, 2ix) 
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such that ip n ,e{x + Z,y) = e t6 ip n) g(x,y). That is, the energy spectrum forms continuos bands 
parametrized by an integer index n and the Bloch parameter 9, and the system eigenfunctions 
are Z-periodic up to a phase 9 = qZ with q the Bloch quasi-momentum. Therefore, we argued 
that wave propagation is ballistic with velocity given by 



V n (6) 



Z dE r! 



h dd' 



(2.27) 



which is proportional to the band slope. Since 



N 



dx5(C-f(x))e[f'(x)}f(x) 



(5.3) 



is the number of points such that f(x) = C with positive derivative in the interval (a, b), Faure 
expressed the number of Bloch modes with energy E propagating with positive velocity as 



N B (E) = J2 d9 5(E- E n (9))Q[v n {9')} 



dE r , 



d9 



(5.4) 



where S(-) is Dirac's delta function and 0(-) is Heaviside step function. Considering the distri- 
bution of quantum velocities Pn{v, E), defined as |AK98j 



P h (v, E) = h 2 J2j^^- <0))5{E - E n {9)) 



(5.5) 



with v n {9) defined in (2.27), one obtains 

N B (E) = 



hL 



vP h (v, E) dv. 



(5.6) 



The normalization of Pn(v, E) over velocities v is 

J dvP h (v,E) = h 2 J ^P(E,9), 



with 



p{E,6) = y £ / 6{E-E n {e)) 



(5.7) 



(5.i 



the density of states with fixed Bloch parameter 9. We note that (5.6) is an exact result, a par- 



ticular case of the Kac-Rice formula [Wil99], valid for quasi-one-dimensional periodic systems. 
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The number of propagating Bloch modes ( 5.6[ ) is a rapidly fluctuating quantity as a function 
of E at scales near the mean level spacing AE = 1/pe, where 



Pe 



d9 
2^ 



(5.9) 



is the mean density of states and the average (-)e is the semiclassical ensemble average defined 
in section 4.1 The fluctuations of N B (E) are eliminated by averaging it over this ensemble. 



Hence, we now consider the average number of propagating Bloch modes, 



(N B (E)) 



1 

hZ 



v(P h (v,E))dv. 



(5.10) 



In what follows, we will drop the average label E unless it is necessary. As follows from (5.7), 
(Ph{v, E)) is normalized to h 2 p B = Ve where v B is the Liouville measure of the constant energy 
surface T, E = {(x,p x ,y,p y ) : < x < Z, H(x,p x ,y,p y ) = E} (Weyl law). The averaged 



velocity distribution in the rhs of (5.10) can be analyzed semiclassically to obtain (Nb(E)}, as 
we show bellow. 

We now turn to obtain the leading order term of the semiclassical expansion of (N B (E)). 
Asch and Knauf [AK98J proved the semiclassical convergence of Pn{v , E) towards the classical 
asymptotic mean velocities distribution P a (v,E), 



P h (v, E) = P a (v, E) when h ->• 0, 



(5.11) 



that holds for test functions independent of h (i.e., over intervals of order Sv ~ O(h ), 5E 



O(h )), a coarse graining consistent with the average taken in (5.10). In the right-hand side of 



(5.11), P a (v,E) is the probability density of the asymptotic mean velocity 



lim ?W. 

t— >oo t 



(5.12) 



with x(t) the longitudinal position at time t of a particle in the waveguide whose initial condition 
was taken randomly with a uniform probability distribution on the surface When the 

classical dynamics is not purely ergodic, so there are some tori associated to ballistic trajectories 
in the stroboscopic Poincare section S x , the leading order term in the semiclassical limit of Ng 



can be obtained just using equivalence (5.11 ) in the integrand of (5.10). In this case, one obtains 



(E) + o{hr 



(5.13) 
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where 



A*bal(-E0 



1 

Z 



v(P h (v,E))dv 



(5.14) 



can be identified |Fau02| as the measure of ballistic trajectories in S x [see equation (3.41)]. On 
the other hand, in case the dynamics is completely ergodic, P a (v, E)dv is a punctual measure 
at v — because, as implied by the CLT (3.29), 



lim 

t— >oo 



1 x(t) 



lim —p 

i^ 00 y/t 



(5.15) 



with x a Gaussian random variable with null mean. Hence, v a — almost surely. In order to 
obtain the leading semiclassical contribution in (Nb(E)), it is necessary to quantify more pre- 
cisely how the quantum level velocity distribution approaches the classical asymptotic velocity 



distribution as K — > 0, that is, we need a refinement of (5.11) 



As discussed in section 3.2.2 when the classical dynamics is fully ergodic and the mixing 
rate is rapid enough, we expect that (almost) any initial ensemble of particles in the system 
will relax by diffusion, with a Gaussian density profile spreading according to (x(t) 2 ) = D B t for 
large times, with De the diffusion coefficient for Eg. Hence, the mean velocity distribution at 
time t is given by 



Pt(v,E) 



exp 



2a\ 



where of = -D^/t. In the limit t — > oo, P t (v,E) converges to a punctual measure at v 



(5.16) 



0. It 



has been shown |EFK+95j . using the Gutzwiller trace formula and also with RMT arguments, 
that in systems with hyperbolic classical dynamics the variance of a quantum operator in the 
semiclassical limit is equal to the variance of the associated classical observable at time t = tjj/g, 
with tjj = h/AE = ve/Ii the Heisenberg time of the unit cell and g a factor depending on the 



anti-unitary symmetries of the system. For the billiard we study in section |3.4| the factor is 
9 

section 



2 due to the transversal mirror symmetry of the unit cell [RB86J [see the discussion in 

(v n (8) 2 ) of the quantum 



5.1.3 



Then, we have a relation between the variance o\ 



level velocity and the variance of the classical asymptotic mean velocity at, namely 



4 



a 



1 H 



for h — > 0. 



(5.17) 



This relation between level velocity variance and the classical diffusion coefficient De was studied 
in detail for the kicked rotor |LCT99] where the diffusion coefficient displays a nontrivial behavior 
as a parameter is varied. In a different context, [BGOlJ explores the relation between classical 
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diffusion coefficient and quantum evolution properties. 



Inspired by (5.17), Faure conjectured a higher order semiclassical equivalence between the 
quantum velocity distribution and the classical mean velocity distribution at time t = t B /g, 



P h (v,E) = P i 



t=t H /g 



'v, E) when h — )■ , 



(5.18) 



over widths of order SE ~ O(h ), 5v ~ 0(\/h). This conjecture is based on result (5.17) and the 



Quantum Ergodicity Theory [BSS98J. It can also be seen as a consequence of the universality 
in the parametric variation of energy levels discovered by Simons and Altshuler [SA93aJ, as we 



will discuss in section 5.1.3 Using (5.18) in the integral of (5.10) we have, for a purely diffusive 
waveguide, 

1 



(Nb(E)) 



Vh\ 2vrZ 2 



(5.19) 



Equation (5.19) was derived in |Fau02| -without the anti- unitary symmetry factor g- and 
checked numerically for the kicked Harper model. 



5.1.2 (Nb) in a periodic diffusive billiard 

Billiards where presented in chapter [3j they are a particular class of Hamiltonian systems 
defined by free point particles which collide elastically with the billiard (waveguide) walls. 
Therefore, they have the simplifying property of equivalent phase-space dynamics on every 
energy shell except for a time scale change. For two-dimensional billiards, v B = 2nA c depends 
only on A c , the area of the unit cell. Hence, the semiclassical limit h — > is equivalent to 
k — > oo due to the relation h 2 k 2 = 2E [see (2.1)]. In the asymptotic mean speed v a takes 
values in the interval [0, y/2E\. Considering the change of variables v = \/2Eu we have from 



(5.10) and (5.11) 



(N B (E)) 



V2E 
~hZ~ 



uP a (u) du 



(5.20) 



with P a {u) the classical asymptotic mean velocity distribution in £1/2 where the particles move 
with constant unit speed v = 1, after setting their mass m — I. Since \/2E = Hk, 



(N B (k)) 



2nZ 



vP a (v) dv. 



(5.21) 



Thus, for ballistic waveguides, i.e. when P a {v) is not a degenerate punctual measure at v = 0, 
the average number of propagating Bloch modes goes asymptotically for k — > 00 as (Nsik)) ~ 
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fibaik/2n + o{k) where \Xhai is defined in (5.14) 



For diffusive billiard chains, the diffusion coefficient satisfies De = \v\D\ with |u| = y~2E 
and Di the diffusion coefficient in the energy shell E1/2. Then, using \v\ = hk and v% = 2nA c 



in (5.19), we obtain that the average number of propagating modes in a diffusive waveguide is 



(N B (k)) = v^S?! y/k + {Vk) for k ->• oo. 
v2iiZ 



(5.22) 



5.1.3 Connection to the universal parametric correlation function 



In section 5.1.1 we closely followed Faure's derivation of the semiclassical expression for (Nb) 



given by (5.19). Here, we provide an alternative argument allowing us to better understand the 



connection of this result with |SA93a[ [DMSS97, DMSS98J. We define the probability distribution 



of quantum level velocities used in section 5.1.1 in natural units w, 

dEr, 



E ) = ^T,[_ deS[E-E n (9)]S 



w 



dO' 



so that 



N B (E) = 2n wp(w, E) dw. 



(5.23) 



(5.24) 



Then, as was noted in [KZ92J, the distribution of level velocities satisfies p(w, E) = lim^o </>K(w<j>, < 
with K(Q, 0) the autocorrelation function 



K(n, <p) = ( P (E + n,e + <t>) P (E, e)) E>6 - p E 



(5.25) 



of the energy density p(E,0) defined in (5.8). We recall that the average is over 6 and over a 



classically small energy interval around E (the semiclassical ensemble). According to Simons 
and Altshuler |SA93aj . 



k{u, x ) = K{n ) <t ) )/p\ 



is a universal function with the scaling 



x = ^C{0)(f) and 
C(0) = pl{{d e E n (e)f). 



(5.26) 



(5.27) 
(5.28) 
(5.29) 
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They showed this result for disordered systems and verified numerically its validity for chaotic 
systems in the semiclassical limit. It follows from their analytical expression for k(u,x) that the 
probability distribution of level velocities p(w,E) is a Gaussian with variance C(0)/p 2 E . Thus, 



from (5.24), we obtain 



(5.30) 



where C(0), defined in (5.29), is system dependent and must be computed, for instance, semi- 
classically. In our case, C(0) is simply related to the average of the squared quantum velocity, 



C(0) 



L 2 



(5.31) 



that is given in (5.17). Then, from (5.30) and (5.31), we obtain the dominant term in Faure's 



expression (5.19) for (Nb(E)). Hence, equation (5.19) can be seen as a consequence of the 



parametric level correlations universalitiy. The validity of this universal correlation function 
for periodic systems was studied by Dittrich et al. [DMSS98J in the semiclassical limit. They 
computed K(Q, <fi) by Fourier transforming their semiclassical expression for the form factor 
identifying a C(0) that agrees with the previous result. 

As a last remark, we note that the universal correlation function k(u, x) is different for GOE 
and GUE systems. For periodic systems, where \ plays the role of the Bloch parameter, time 
reversal symmetry T is broken for all x besides three exceptional points [DMSS97J. Thus, when 
no further anti-unitary symmetries exist, k(u, x) moves between GOE and GUE at different 
parts of the band. For the periodic cosine billiard we consider in our numerical calculations, the 
unit cell is invariant under the operator S x that transforms x — > — x (mirror reflection along the 
transversal direction). So, even though the unit-cell Hamiltonian do not commute neither with 
S x nor T, it does with the operator TS X which is anti-unitary. Then, as x is varied, the system 
always belongs to GOE [RB86J . The presence of this anti-unitary symmetry explains the factor 



g — 2 in (5.19) and (5.22), which enhance (N B ) by a factor \/2. 



5.1.4 Numerical results in a periodic cosine waveguide 

We now study numerically the semiclassical behavior of (N B (k)) in the cosine billiard chain. 
In order to compute the transmission and reflection matrices we employ the admittance multi- 
modal method described in |PaglO| . For each value of k, the field is described using the local 
waveguide transverse modes truncated to a finite number of channels, as described in chapter 
M The number of propagating modes N B {k) is obtained by counting the solutions of the unit 
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cell transfer matrix eigenvalue problem 



M 1 ^ = \ip, (2.25) 



with |A| = 1 and positive group velocity. As we have discussed in this chapter, N B (k) fluctuates 
very fast as a function of k at scales Ak = 27r/(kA c ) [] We consider an average of the form 

(N B (k)) r = J N B {k')f r (k - k') dk' , (5.32) 

where f r (k) is a positive function normalized to one with compact support in an interval of 
length rAk with r ^> 1. The dependence of (NB{k)) r in the width of f r (k) is discussed at the 
end of this section but for what follows we drop explicit reference to r. 



In figure 5.1, we show (N B (k)) for k G [tt, 1507r], for a cosine billiard with A\ — 1 and 



A2 = 0.3 [see section 3.4 for the cosine billiard definition]. As we have discussed in section 3.4 
for this billiard configuration the phase-space does not exhibit noticeable tori and its dynamic 
is strongly diffusive. In the k interval explored, the number of open transverse modes in plane 
leads ranges from 1 up to 150. For all wavenumbers, 100 evanescent modes where used in 
the numerical calculation. The dashed line shows a Gaussian smoothed moving average of 
N B {k) with variance tt 2 for a uniform sampling of wavenumbers with spacing 5k ~ 0.377T, 
which is bigger than the mean level spacing in the whole k range. The full line is the expected 



semiclassical result (5.22) adjusted with an additive constant. This constant appears because 
in the system under study (N B (k*)) = for some k* > and we are not sufficiently deep in 
the semiclassical regime for its relative value to be negligible. N B (k) fluctuates very rapidly 
with k, as expected, but the agreement between the moving average and the semiclassically 
predicted curve is quite good in the whole wavenumber interval explored. We observe similar 



good results for all other configurations of the unit cell discussed in section 3.4, even when 
there are some small stable tori in the S y Poincare section. This is not surprising because the 
classical dynamics in those systems is also strongly mixing, and for waves to resolve such small 
structures wavenumbers of order k > 20007T would be needed, a fairly large value compared 



to the range considered. In the inset of Fig. 5.1, we show the diffusion coefficient calculated 



classically (from Montecarlo simulations) and also the one obtained from the best fit following 



(5.22) for the quantum computation of N B {k) for five values of A 2 . 



As we mentioned in the definition of the average number of propagating Bloch modes (5.32 ) 



the width of the interval rAk over which we take the average is arbitrary and it is only assumed 



'Ak follows from the Weyl law N = k 2 A c /(4w) 
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Figure 5.1: The dashed line is the Gaussian smoothed moving average (with standard deviation 
~ 7r) of the number of propagating Bloch modes Ns(k) in the cosine billiard chain with Ai = 
1 and A 2 = 0.3 for a uniform k sampling with increments 5k = 0.371477T. The expected 
semiclassical result for (A^^(fc)) (full line) shows good agreement with the moving average. The 
four black points (with their associated standard error) were computed for the same geometry 
using a uniform sampling of 153 wavenumbers in an interval 300 mean level spacings wide [see 
figures 5.2 and 5.3) , which constitutes better statistics than the moving average (a sampling at 
least 15 times more dense). In the inset, the diffusion coefficient calculated classically (black 
dots) and the one obtained from the best fit following (5.22) for the quantum computation of 
Ns{k) (circles) for five values of A 2 is shown. In all cases, a good agreement with the expected 
semiclassical behaviour given by (5.22) is observed. 
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Figure 5.2: (NB{k)) r as a function of r for four energies in cosine billiard with A\ = 1 and 
A2 = 0.3. A uniform sampling of 153 wavenumbers distributed in a 300 mean-level-spacing 
wide interval was used. In all cases after r ~ 200 the average fluctuations are notably reduced 
and the result for (Nsik)) is in agreement with the expected semiclassical result [see figure |5\T] 
indicating that this is a good smoothing window to eliminate the spectrum fluctuations. 

to be much larger than the mean level spacing Ak and classically small. In order to find out 
how relevant this choice is, we calculate (NB(k)) r for four different values of k as a function 
of r G [2, 300] using a uniform k sampling of 153 points around the four wavenumbers. In 
figure |5.2 we see that, after a transient phase, at r ~ 200 the value of the average reaches an 



approximately stable region in the four cases. The four black dots in figure 5.1 show [Nsik)) 



for r = 300 and the k values of figure 5.2, displaying better agreement to the theoretical curve 
since this average takes more semiclassical realizations. 

We end this section remarking how deep into the semiclassical limit the presented results 
are. In our calculations A//lo -the wavelength in units of the waveguide width- is as small as 
0.012. We note that the calculations in [DMSS97J, which are also compared with semiclassical 
results for periodic billiards, were computed with the parameter X/Hq ~ 0.12. 
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Figure 5.3: Number of propagating Bloch modes Nb as a function of k for the cosine billiard 



with Ai = 1 and A2 = 0.3 for the four fat points shown in figure 5^ Each fat point is the 
average of the 153 uniformly distributed wavenumbers k shown in the four plots. 

5.2 Random matrix periodic chain model 

5.2.1 (Nb) for a random matrix periodic chain 

A periodic waveguide random matrix model was introduced in chapter |4j It consisted of a 
periodic chain of strongly scattering cavities whose S matrix is taken from the COE or CUE 
ensembles, depending on the desired symmetry of the system. In this section, we obtain an 



expression for (Nb) in such systems by means of Fame's result (5.22) and the Machta-Zwanzig 



diffusion coefficient approximation (3.37). 



In chapter [4j we discussed how the circular ensembles, COE and CUE, characterize the 
statistical properties observed in a quantum dot scattering matrix. We defined a quantum dot 
as a coherent strongly scatterering cavity, with chaotic classical dynamics and such that its 
mean residence time r is much larger than its ergodic time. From these properties follow that 



the S matrix is insensitive to the microscopic details of the quantum dot [see equation (4.24) 
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Figure 5.4: (Right) Number of propagating Bloch modes (Nb) calculated numerically for 
the COE (circles) and CUE (squares) periodic chain model; the black lines show the the- 
oretical result given by (5.36). (Left) The circles show the log 10 of the relative error e = 



(N B ) 



(numeric) 



/(Nb) 



(thooric) I 



as a function of log 10 (iV) for both ensembles and the line, which 



have a slope —1, indicates that a = 0.5 [see equation (5.36) 



and hence implies universality in the sense of RMT. Not by coincidence, as we argued in section 



4.3.3, the conditions under which the Machta-Zwanzig diffusion coefficient approximation (3.37) 



holds are the same and thence it gives the universal diffusion coefficient for (the ray dynamics) 
in a chain of quantum dots, which depends only on the unit cell perimeter and openings length. 



Let the quantum dot periodic chain be defined as in section 4.3.4, with cavities of area A 



and opening of width ho. Then, its Machta-Zwanzig diffusion coefficient (3.37) is given by 



D 



MZ 



2Z 2 h 
7rA r 



(5.33) 



We know that for a periodic chain of chaotic cavities in the semiclassical limit, the average 



number of propagating Bloch modes (Nb) is given by (5.22). Hence, inserting (5.33) into this 
expression we obtain 

(Nb) = (5.34) 



7T- 



Recalling that the number of open modes in the leads is given by 



N= [h k/7r\, 



(2.5) 
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we finally have the average number of propagating Bloch modes in a periodic quantum dot 
chain, 



(Nb) 



\- + o(Vn) 

V 7T 



for N — > oo. 



(5.35) 



In figure [5T4| we compare this result with montecarlo calculations for the COE and CUE models 
showing perfect agreement in the CUE case and a constant additive correction —0.2 for the COE 
case. The latter can be attributed to a weak-localization effect, as we will argue in chapter [6} 



We also studied the next order correction to (5.35), which turned out to vanish as ~ iV l l 2 for 



large N for COE and CUE as is also shown in figure 5.4 Therefore, we have that 



- 0.2 + 0(N- 1/2 ) for COE 



(Nb) = { 



(5.36) 



+ 0(N' 



-l/2s 



for CUE. 



5.2.2 Characterization of Nb distribution 

As we have seen in the previous sections, once the semiclassical ensemble or RMT model is 
introduced, Nb can be considered as a random variable taking integer values for each realization 
of the waveguide. We now study numerically some aspects of its probability density function, 



n 



using the RMT periodic chain ensemble. 



In figure 5.7, we show the pdf histograms for COE and CUE with N = 10,50,100 and 



compare them with the best fitted binomial distribution imposing that (Ng) follows (5.36); the 
agreement is good in both cases even for rather small t matrix dimension N = 10. The binomial 
distribution is the discrete pdf of the number of successes in n independent realizations of a 
dichotomic experiment with success probability p; its mean and variance are given by np and 
np(l —p), respectively. In our case, the analogy works taking n = N, p = (Nb)/N = 1/y/nN 
and each success meaning the presence of a propagating Bloch mode in the spectrum. However, 
this is not a very good approximation because the experiments are not independent. In fact, 
given there is a Bloch mode in the spectrum, the probability of a second one is not the same as 
for the first because of correlations. As we see below, this fact causes the variance of Nb to be 
less than the expected from the binomial distribution, or equivalently, we have to set n ^ N in 



the fitted binomial distributions shown in figure 5.7 
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Figure 5.5: Variance of the number of propagating Bloch modes Nb as a function of the number 
of open modes in the leads N, in a RMT periodic chain model with the COE (circles) and CUE 



(squares) ensembles. The gray lines show the fit given by (5.37). 



The variance of Nb is shown in figure 5.5 and is found to follow 



Var[iV 



B\ 



N 

0.82W— -0.24 



0.41 



for COE 



for CUE. 



(5.37) 



We note that this means (Nb) /Nsx[Nb\ — > a with a = 0.82 in the COE case and a = 0.41 in 
the CUE case, for N — > oo. The fact that (Nb) /Vsix[Nb\ converge to a constant value remind 
us of shot noise [dB96j with Bloch modes playing the role of the quantized packet (electron, 



photon, etc). Since conductance g ~ Nb as we discussed in section 2.2.2, the conductance in a 



long chain will display this shot noise fluctuations [see section 6.4 



To finish this section, we consider P^ B (0) as a function of N. This quantity gives us a 



measure of the average band gaps density in the ensemble. In figure 5.6 we show that P/v s (0) 
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Figure 5.6: Logarithm of the probability of having zero Bloch modes Pjv s (0) as a function of N 
in the RMT periodic chain model for COE (circles) and CUE (squares). The gray lines show 



the fit given by (5.38). 



is well described by 




Pn b (U)= { 



exp 



exp 0.37-1.8 




for COE 



for CUE. 



(5.38) 



Surprisingly, expression (5.38) for the COE case agrees exactly with what would be expected 



from a Poisson distribution, in which the probability of observing zero events (in our case zero 
Bloch modes) is e~ A with A the distribution average. The Poisson distribution is a binomial 
distribution in the limit n — > oo and p — > 0, which is consistent with our case where n ~ N and 
p ~ 1/ y/N. Of course, again the analogy is broken for the rest of the distribution because of 
the non-null correlations in Bloch spectra. 
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Figure 5.7: iV^ Probability density function (bars) for a COE and CUE periodic waveguide 
with iV = 10, 50, 100. It is also shown the best fitted binomial distribution (dots) with mean 



(Nb) [see equation (5.36)] and variance given by (5.37). 
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Chapter 6 

Conductance in diffusive periodic 

waveguides 



In this chapter, we study the conductance of periodic quasi-one-dimensional waveguides with 
diffusive classical dynamics. The scattering approach to waveguides was presented in chapter 
[2j where we defined the Landauer conductance which is the main focus in our analysis. On the 
other hand, we have discussed relevant classical properties of chaotic billiards in chapter |3j and 
introduced a random matrix model of periodic diffusive waveguides in chapter |4j We employ 
this model to obtain statistical properties that hold in general diffusive periodic waveguides. 

We consider a periodic waveguide with L unit cells connected to two semi-infinite plane 



leads at its extremes as defined in section 2.1.1 The classical dynamics of unit cell is assumed 



diffusive, i.e. chaotic, with finite-horizon and strongly mixing. We will take averages over the 
two ensembles defined in chapter |4j the first is given by the RMT periodic chain model which 
we denote (•) (we use the same notation for the ensemble average in disordered systems but the 
context will distinguish which we are referring to), and the second by the semiclassical ensemble 
which we write (•}*.. 

Classically, for disordered or periodic chaotic systems, Ohm's law states that the conductance 
scales as g ~ L d ~ 2 where d is the dimension, so for a wire (d — 1), we have g ~ 1/L. In quantum 
systems, the conductance also scales in this way if wave coherence is broken at lengths of the 
order of the mean free path. Otherwise, in a wire with disorder of any strength, the average 
conductance decay exponentially (g) ~ e~ 2L ^ for long wires, due to the Anderson Localization 
effect [ATAF80J. This is also observed for d — 3 in the presence of strong disorder but in 
two-dimensions the decay is not necessarily exponential. The conductance exponential decay is 
triggered for wires longer than the localization length £ = lim^oo — 2(<9(log g(L)) / dL)~ l and for 
shorter wires an Ohmnic regime (g) ~ 1/L is still present, which is known as the metallic regime 
and characterized by the presence of weak localization and universal conductance fluctuations 
|LS85] . 

In contrast, it is clear that in quantum periodic waveguides there is no localization (outside 



band gaps) because of Bloch ballistic modes. In section 6.1, we study the conductance of long 
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Figure 6.1: Conditional mean (An b \Nb) for six COE samples of 100000 realizations for dimen- 
sion N — 25, 50, 75, 100, 125, 150. 



periodic chains and establish in more detail their dependence on Nr as proposed in section 



2.2.2 Then, in section 6.1.2 we show that the non-propagating part of the conductance has an 
asymptotic long-chain limit that decays as L~ 2 . The main result of this chapter is the existence 
of an ohmic regime in periodic diffusive waveguides as discussed in section |6~2| Finally, we study 



the conductance fluctuations and distribution in the ohmic and Bloch-ballistic regimes. 



6.1 Long-chain conductance properties 



6.1.1 Number of Bloch modes and conductance average 



As seen in section 2.2.3, in a chaotic waveguide with periodic structure the typical conduc- 



tance, i.e. the conductance for a fixed energy taken from the semiclassical ensemble, displays 
exponential decay towards asymptotic quasi-periodic oscillations in the long chain limit [see 



equations (2.38) and (2.46)]. From (2.34), we can write the conductance as 



g(L)= 9oD (L)+ ]T Tt{L), 



(6.1] 
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where 



9oo (L) = T ^ L )- 

|Ai|=l 



(6.2) 



Note that the last term in the rhs of (6.1) decays exponentially [see section (2.2.2)], so g(L) 
(7oo(L) for L ^> £, with 



(2.37) 



min{log|Ai|} 

|A 1 |>1 

the largest evanescent mode decay length, defined in chapter [2j The asymptotic long-chain 



conductance g^L) is quasi-periodic in L as can be inferred from (2.46). Also, note that all Tj 
and hence g depend on the ensemble realization, this is, on the wavenumber k or on the RMT 
realization. Then, taking goo{L) semiclassical average one obtains a constant value, 



(g<x>(L -> oo)) fc = (goo) k < (N B )i 



(6.3) 



as is plotted in figure 6.3 In order to make the inequality in (6.3) more precise, we write the 

goo(L) = N b A Nb 



long-chain conductance as 



(6.4) 



with 



A 



N B 



l N B 

-^T t (L) forA^O 



N B ~ 1 



(6.5) 



1 for N B = 

the average Bloch transmission, which is a random variable (of the semiclassical or RMT en- 



semble) dependent on N B (and L) and bounded between in [0,1]. Figure 6.1 plots how the 
conditional mean (An b \N b ) decreases as a function of iVe from a value close to 1 for N B = 1 
and saturates around 1/2 for large N B . Hence, we have 



N 1 
(goo(L)) = ^2(A Nb \N b = n)nP NB {n) > -(N B ), 



(6.6) 



n=l 



where P^ B {n) is the probability of having n Bloch modes [see section 5.2.2 . Therefore, from 



(2.38) and (6.6) follows that 



(N B ) > (g(L)} > -(N B ) forL^oc 



(6.7) 
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in the semiclassical limit. Numerically, it is observed using our RMT model that a good ap- 
proximation for large N is given by 

hm( 9 (L)) = / 07 » OE - a20 f ° rC0E (6.8) 
[ 0.75 (iV B )cuE- 0.15 for CUE. 

For the cosine billiard this result is also good and in fact should be universal since it stems from 
the RMT Ti eigenvalues distribution. The difference in the conductance for COE and CUE can 



be seen as consequence of weak localization as discussed in section 6.3 



6.1.2 Approach to the asymptotic regime 

The decay rate of the conductance g(L) to its asymptotic regime depends on I [see equations 



(2.38) and (2.37)]. We begin this section focusing on the conductance non-propagating part 



g np (L) = g(L) - 9oo (L) = ]T T,(L), (6.9) 
which is dominated by the slowest to decay evanescent transmission amplitude, 

r " (L) = 2 + A„ W 4 + A mW - ^ 4am(LrV2t/ ' f ° VLZe - ( " 0) 



where a m (L) _1 ~ 0(1) is quasi-periodic as we discussed in section 2.2.2. Hence, taking the 
average we can drop a m (L) since the decay will be dominated by the exponential, this is, we 
have 

(g np (L)) = <4e- 2 ^>- (6.11) 

In order to calculate this average we will study the pdf of the decay length £. We can inspect nu- 
merically this distribution using RMT; we focus on COE since this is the ensemble representing 
the cosine billiard without a magnetic field. We consider the pdf 

P(£,N)= [ dfi CO E(U T U) 5 [I - (max|£(C/ T C/)|)] , (6.12) 

JUeU(2N) 
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0.001 



Figure 6.2: Scaled pdf P(£) of the longest decay length i for the COE ensemble (gray line), 
constructed by superposing the pdf N 3 I A P(W 3 ' A ,N) for N = 25,50,75,100,125,150. The 
dashed black line decays as £~ 3 . 

where /zcoe is the invariant COE measurd^ and 



C(S) = {(log (A))- 1 : A e eig(M(5)) A |A| > 1}, 



(6.13) 



is the decay lengths set associated to the spectrum of the transfer matrix M(S) obtained from 



a scattering matrix S. We have found an asymptotic scaling, as shown in figure 6.2, given by 



N 3/4 P(£N 3/ \ N) — ► PU) 



(6.14) 



with the scaled pdf P(£) displaying an algebraic tail, 



P(t) ~ I' 3 for i -»■ oo. 



(6.15) 



Horvat and Prosen[HP07] studied a pdf similar to (6.12) but for the maximum modulus M 
eigenvalue and found and analogous asymptotic scaling with exponent 1/2. Finally, using the 



*The 2N-dimensional COE ensemble measure /icoE is the induced probability measure in the symmetric 
space U(2N)/0{2N) by U(2N) Haar measure. 
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Figure 6.3: Average conductance (g(L))k (squares) and its non-propagating part (g n p{L))k (cir- 
cles) for the cosine billiard chain with A\ = 0.5, A 2 = 4.5 and k = 30.337T. The average number 
of propagating modes {NB)k = 2.6 (dashed line) is an upper bound for (g(L))k as L — > oo [see 
equation (6.7)]. On the other hand, (g np (L))k displ ays an algebraic decay ~ L~ 2 (gray line) for 



large L as predicted by the pdf P(t) [see equation (6.16) 



scaled pdf (6.14) and its algebraic tail (6.15), we have 



(g np (L)) = <4e- 2 ^> ~ IT 2 . 



(6.16) 



This result, which was obtained using the RMT periodic chain model, also holds for the semi- 



classical ensemble as expected and can be seen in figure 6.3 for the cosine billiard chain. 



In disordered systems, a standard characterization of the conductance decay to the asymp- 
totic regime (which in that case is localization with (g) — > 0) is the linear decay of the av- 
erage logarithmic conductance (logg(L)) ~ — L for long wires. Interestingly, although in pe- 
riodic chains the asymptotic conductance is not null because of the presence of propagating 
Bloch modes, in diffusive chains (log g(L))k shows a similar behavior. In fact, we have for 
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<log(g)> 
2.0 




-1.0 



Figure 6.4: Average logarithmic conductance (log g(L))k (dots) for the cosine billiard chain with 
A\ = 0.5, A 2 = 4.5 and k = 30.337T as a function of L. The linear decay ~ — predicted by 



(6.18) is clearly observed (gray line). In this case £ = 98.6, P/v s (0) = 0.06 and N = 15 hence 
(6.20) implies the harmonic mean decay length (i^ 1 )^ = 0.83 (in unit cell length units). 



L » ( 9oo ) k 1/2 , 



(log g(L)) k = P Nb (0) (log4e- 2L ^) fc + (1 - P Nb (0)) (log 9oo (L)) k 
w -(2^ 1 ) fc P Nb (0) L + constant 



(6.17) 
(6.18) 



assuming correlations between and i are negligible. Therefore, even though in a periodic 
chain there is no exponential decay of the average conductance and its limit is non-null, 



(log g{L)) k 



2L 

T 



for large L 



(6.19) 



just like in a disordered wire. However, the constant which corresponds to the localization 
length in disordered wires, in periodic chains can be seen as a measure of the band-gaps in the 
spectrum, since it is proportional to the probability of having zero Bloch modes, 



(6.20) 



where (-^ is the average taken over the scaled pdf P(£). The probability Pjv s (0) as a function 
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of N was studied in section 5.2.2 for the COE and CUE ensembles. We have numerically 
estimated the scaled harmonic average decay length to be = 0.84 using the scaled pdf 



(6.14). Alternatively, we can compute this number indirectly measuring £ 1 and using (6.20). 



In figure 6.4, we plot (log5' 00 (L)}fe in the cosine billiard chain and use £ and P/v s (0) to calculate 



the scaled harmonic average decay length 
l/g(L) 



/ oo 




Figure 6.5: Average resistance (l/g(L))k (black dots) for the cosine billiard with A\ = 0.5, 
Ai = 4.5 and k = 30.337T. The full line shows the ohmic regime holding for L < 10. Also, 
several realizations of l/g(L) are displayed in gray: a few cases with Nb = grow exponentially 
and the rest reach an asymptotic oscillating value l/p co (L) ~ 0(1/Nb)- The exponential grow 
of the average resistance (R)k for L > 10 is explained by the presence of the rare cases with 
Nb = 0. 
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<N/g) C0 E 




L 



Figure 6.6: Plot of the average normalized conductance (g(L)/N) and resistance (N/g(L)) 
(inset) for the COE periodic chain model with N = 10, 20, 30, 40, 50 (dots, squares, triangles up, 
triangles down, diamonds) as a function of L. The average conductance collapse to (g(L)/N) ~ 
(L + l)" 1 and the average resistance to (N/g(L)) ~ L + 1 (full lines) for large N during the 
ohmic regime which is observed in a range 1 < L < \/N [see discussion around 

6.2 Ohmic regime 

As we have seen, in periodic waveguides, even if its classical dynamics is chaotic, quantum 
transport is ballistic due to the existence of propagating Bloch modes. Interestingly, however, 
this asymptotic ballistic regime is preceded by an ohmic regime for short chains, similar to 
what is observed in disordered system. As we have discussed in chapter [TJ an ohmic conductor 
shows the conductance scaling G = cxW/L, where W and L are the sample width and length 
respectively, and o is its conductivity. In mesoscopic systems, the coherent propagation of 
waves causes this scaling to break up at large lengths; in disordered systems there is an ohmic 
regime for L <C £ followed by localization and, as have been discussed in this work, in periodic 
systems propagation is ballistic in long-chains with g(L — )■ oo) ~ A?£. In this section, we study 
the existence of the ohmic regime in diffusive periodic waveguides. We consider the resistance 
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defined as 



R= m (6.21) 



which is plotted in figure 6.5 for a cosine waveguide. Qualitatively, we found Ohm's law (R) = 
(L + 1)/N holds for small L, i.e. before (g(L)) reaches its asymptotic value. This regime is 
followed by localization-like exponential grow of the resistance. The latter may seem surprising 
but is a consequence of the non-null probability of = (or equivalently P^ B (0) ^ 0). In 
fact, although Pjv b (0) ^ decays to zero as iV — > oo, it dominates in the average (l/g(L)) 
for long chains. Of course, for particular realizations of l/g(L) the most probable is l/g(L) = 
l/#oo(£) < oo with N B ^ 0. 

On the other hand, the existence of the ohmic regime is a consequence of the full M matrix 
spectrum and is observed for small L, where several non-propagating Tj still contribute to the 



sum in the rhs of (6.1). In order to understand this result we digress briefly to disordered 
systems, where the linear scaling of R(L) observed before localization sets in can be explained 
by the DMPK localization lengths spectrum £ n , n — 0, . . . , N — 1. It is known [Dor 8 2] that 
in(N)^ 1 ~ n/N£, where N£ denotes the disordered wire localization length. From the polar 



decomposition (2.31), the disordered wire conductance is given by 

N ~ 1 4 

g = V s ^. (6.22) 

^ 2 + e L /^ + e-V^ 

Then, the ohmic scaling 1/g ~ L for £ <C L £N in disordered system can be seen as a 
consequence of l~ x lineal dependence on n. It turns out that the spectrum |A n | of Mi = Mi(S) 
with S taken from COE has a similar property, namely f" 1 = log |A n | ~ n/N for (iVs) < n ^ N 
which explains the existence of an ohmic regime in diffusive periodic waveguides [see figure 6.7 . 
This is a statistical characteristic of the transfer matrix spectrum, which can be recasted stating 
that 

P(|A|)~pq-, for ^<log|A|«l (6.23) 

withp(|A|) the marginal pdf of the spectrum absolute values |A n |. From this follows that 1/g ~ L 
for L < N/ (Nb) ~ \/N, which can also be understood from the two following arguments. 

First, one expects the ohmic regime to be observed only for chains whose diffusion time 
to = L 2 /D is smaller than the unit cell Heisenberg time tu = v^jh at which Bloch states 
are resolved, where we recall that D = Div, with v the particle's speed, and for a billiard 
Ue = khA c /v. Then, we obtain that to < tn if and only if L < (Nb) ~ VN which is the 
expected result. Note that (N B ) = aN/(N B ) with a = A c D 1 /(h Z 2 ) ~ 0(1) so there is no 
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contradiction with our original estimate. The second argument comes from the matching of the 
diffusive regime conductance (g(L)) & g /(l + L), with go = N, to its asymptotic average value 
(9oo) ~ (Nb) from where we obtain again that the ohmic regime is confined to L < N/(Nr)- 

The ohmic regime of diffusive periodic systems is of course also observed in our RMT periodic 
chain model as plotted in figure 6.6[ where we see that 



(g(L)/N) = (N/g(L))-i = (L + l)-\ 



(6.24) 



We note that this is one of the signatures of the metallic regime in disordered wires and a 



consequence of the 0(1) fluctuations of the conductance. We come back to this in section 6.5 



Note that we have shown that Ohm's law in periodic diffusive waveguides is a property of the 
typical resistance R = 1/g (without ensemble average) in contrast to its exponential regime 



which is observed in a set with measure going to zero as N — > oo [see figure 6.5 




Figure 6.7: Average inverse decay lengths spectrum l~ l = log|A„| for the COE periodic chain 
model. The plot shows Ni~ l for N = 10, 20, 30, 40, 50, 60, 70 (circles, squares, up triangles, down 
triangles, stars, diamonds). The integer n is the absolute- value-sorted A n eigenvalue index. For 
this system (Nb) = 4.7. As can be seen in the plot Ni' 1 is close to zero for n < (Nb) and is 
followed by a range with linear growth Ni" 1 ~ n (dashed line) which explains the existence of 
the ohmic regime in the periodic chain for L < (Nb)/N. In the inset, part of the pdf P{N[^ 1 ) 
for iV = 70 is plotted, showing a constant range for 10 < Ni' 1 < 30 which is equivalent to 
PflAl) ~ | | — 1 [see rijjjl]. 
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Figure 6.8: Plots of <5g = (g(L))coE — (<?(£)) cue for AT = 10, 20, 30, 40, 50 (dots, squares, tringles 
up, tringles down, diamonds) showing the weak localization correction in a RMT periodic chain 
as a function of its length L. 



6.3 Weak localization 



As we discussed in section 4.3.3 , the average conductance in a disordered wire or quantum 
dot (a single cell chain) differs for systems with and without time reversal symmetry, namely 
1>9 = (g)coE - (9}cve is given by 



4 
1 



for a quantum dot 



for a disordered wire. 



(6.25) 



This difference, called the weak localization correction (WLC), is a purely quantum effect which 
stems from the destructive interference of time-reversed trajectories present in time-reversal 
symmetric systems (corresponding to the COE). 

On the other hand, as we saw in section 5.2. 1[ the number of propagating Bloch modes in 
the RMT periodic chain model also displays a weak-localization-like effect shown by equation 
(5.36), this is SN B = (N b ) C oe — (-Ab)cue = —0.2. Since in long periodic chains (g) ~ (Nb), 
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then we also expect a WLC in the Bloch-ballistic regime. This is shown in figure 6. 8| where Sg 



is plotted as a function of L for the RMT periodic chain model. By construction, we expected 
to recover Sg(L = 1) = — 1/4 (the quantum dot WLC) as shown. Then, for 1 < L < y/N, i.e. 
during the ohmic regime, 5g undergoes a transient stage before approaching a constant plateau 
close to —0.2 for L ^> \/N, i.e. in the Bloch-ballistic regime, as expected. We note that the 
magnitude of Sg is largest during the ohmic (metallic) regime taking a value ~ —0.3, which is 
close to the WLC in a disordered wire. 

In conclusion, although the WLC is usually associated with the metallic regime, here we 
see that a small difference 8g = —0.2 also persists in the Bloch-ballistic regime for long peri- 
odic chains. Recently, the WLC for a periodic system was studied |Tia09] assuming that the 
Ehrenfest time is larger than the ergodic time. In our case, where RMT is a good model for the 
periodic waveguide, the Ehrenfest time is smaller than the ergodic time. 

6.4 Conductance fluctuations 

Conductance fluctuations in chaotic periodic chains, characterized by the variance Var[g], 
behave quite differently to the disordered wire. For the latter, in the metallic (ohmic) regime 
Var[g] ~ independent of N, L and disorder properties, which is the reason they are 

called universal conductance fluctuations. In a quantum dot Var[g] = 1/8 whereas for a 
disordered wire Var[g] = 2/15 /9 _1 , the discrepancy being a result of non-geometric correlations 
of the eigenvalues T n in the DMPK equation[Bee97j. For a long wire, once the localized regime 
is reached the conductance fluctuations diverge, growing such that Var[logg] ~ L. 



Alternatively, as depicted in figure 6.9, the conductance variance in a chaotic periodic chain 
increases linearly with L during the metallic regime and is well described in the range 2 < L < 
0(VN) by 

r , N1 f 0.05 L + 0.09 for COE 
Vax[g(L)] = \ 6.26 
| 0.03 L + 0.04 for CUE. 

By construction, for L — 1 the quantum dot value is obtained. The linear regime given by 



(6.26) breaks at L ~ 0{y N) signaling the diffusive to Bloch-ballistic transition and is followed 



by Var[g(L)] reaching an asymptotic constant value of the same order. This asymptotic value 



can be understood using decomposition (6.4) and the law of total variance to obtain 



Vax\g(L)) = (N b V&t[A Nb \N b )) +Vht[N b (A Nb \N b )} 

< ^(N B ) +Var[N B ], (6.27) 
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where we used that Vsh[An b \Nb] < 1/4 since An b has its support in [0, l£J We know from 
chapter |H| that (N B ) ~ Var[N B ] ~ VN for large N [see fl5.22| ), ( |5.36| ), ( |5.37| )]. As a result we 
have that 

lim Var[#(L)] ~ v 7 ^. (6.28) 



This is shown in figure 6.10 for the COE and CUE periodic chain models. 



^In general, the variance of a random variable with compact support on [a, b] is bounded above by (b— a) 2 /4. 
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Figure 6.9: Conductance variance for the COE and CUE periodic chains as a function of L. 
The p lots show ensembles with iV = 1 to 50. For L < yiV the variance grows linear with L [see 
(6.26)] and reaches an asymptotic value ~ \fN [see (6.28) . 
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Figure 6.10: Log-log plot of the conductance variance for the COE and CUE periodic chains as 
a function of N. The plots show the variance for different values of L (color lines). For large 
N and L > \/N the variance reaches the Bloch-ballistic regime with Var[g] ~ (black line) 



[see equation (6.28)] 
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6.5 Conductance distribution 

In this section we summarize some qualitative properties of the conductance pdf P(g) in 
periodic chains. For short systems in the metallic regime, before the asymptotic Bloch-ballistic 
state is reached at L ~ \/N, P(g) is Gaussian -like the metallic regime of a disordered wire- 



with mean and variance given by (6.24) and (6.26), respectively. Although the variance grows 



with L, there is a range where it is small enough so that (1/g) ~ l/(<?)- 111 fact, consider the 
expansion 

(g- 1 ) = (g)' 1 - VarM <£>- 3 + O {(g)- A ) . (6.29) 
Since (g) = N/(L + 1) we have that deep in the ohmic regime, where L <C \/N < N, Var[g] <C L 



and Var[g] (g)~ 3 <C (g) 1 - Hence, keeping the dominant term in (6.29) we obtain (1/g) ~ 
l/(g)- Note that at the diffusive to Bloch-ballistic threshold L ~ y/N, Var[g] ~ \fN therefore 
apparently the first term is still dominating becasue Var[g] (g}~ 3 ~ 1/N < (g)^ 1 ~ 1/yN. 



However, the approximation breaks down because P(g) losses its Gaussian shape so the higher 



order moments in (6.29) cannot be neglected. 

In figure 6.11 the distribution P(g) of a COE periodic chain is plotted for different values of 
L. P(g) starts being a Gaussian with its center (g) = N/(L + 1) decreasing with L until (g) ~ 
(Nb) ~ VN at L ~ \/N, where the pdf begins to spread and develop a multimodal structure 
signaling the diffusive to Bloch-ballistic transition. The asymptotic long-chain distribution 
converges to this shape with multiple peaks which arise as a consequence of the Bloch modes 
ballistic propagation. These peaks correspond to realizations of the conductance g = NbAn b 
with different iV^- As we have discussed, (An b \Nb) is close to 1 for small Nb and tends to 1/2 
for larger Nb, hence the better defined peaks at lower Nb- This is similar to what is observed 
in CUE chains. 
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Periodic diffusive waveguides 
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Figure 6.11: Conductance pdf of a COE periodic chain model with N = 50 for several values 
of L. We see that up to L = 9 the distribution is clearly gaussian (black line) and at L = 17 
peaks start to develop. At L = 50 the distribution has reached its long-chain stationary form 
consisting in a multimodal structure which arises from the asymptotic quasi-periodic behavior of 
g around its average ~ Nb- Given an integer n the peak just below g = n is due to realizations 
with Nb = n. Note that for Nb = 1 the peak is very narrow and its center close to one but for 
bigger Nb the peaks start to spread and locate further away from its respective integer upper 
boundary. This can be understood from the behavior of (An b \Nb) as a function of Nb [see 



(6.4) and figure 6.1 . 
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Chapter 7 

Conclusions 



This thesis has dealt with semiclassical transport in periodic chains of chaotic cavities. Our 
analysis has been based on the study of Landauer conductance as a function of the chain length 
and energy. For this purpose, we have employed a realistic quasi-one-dimensional waveguide 
model were we performed direct numerical calculations as well as an RMT periodic chain model. 
The link of RMT to chaotic systems is well known[BGS84j and has been mathematically proved 
for closed systems[MH B + 04"t |HMA + 07 |; this connection explains the famous universality ob- 
served in the quantum dynamics of chaotic systems as predicted by the Gaussian ensembles of 
Hamiltonians. The application of RMT to open chaotic systems is more recent|BS90, JPB94J 
and known to be appropriate in cases with strongly chaotic unit cells where the Ehrenfest time is 
much smaller than the ergodic time. In the opposite case, non-universal behaviour is expected. 
We always assumed that systems were in the RMT regime in this work. The connection be- 
tween RMT and semiclassical averages is given by taking an ensemble of energy realizations 
over several mean level spacings, which is a large interval for quantum mechanics but classically 
infinitesimal. We defined our RMT periodic chain model by taking scattering matrices from 
one of the circular ensembles (depending on the symmetries of the system being modeled) and 
composing it L times to obtain a finite periodic chain. 

The system under study consisted in a quasi-one-dimensional waveguide with a scattering 
region made of a chain of L identical cavities connected to two semi-infinite plane leads at its 
extremes. We assumed the underlying wave equation was the Helmholtz equation, i.e. waves 
(and their corresponding classical particles) were free inside the waveguide and only scattered 
by the geometry of the boundaries. We focused on cavities whose classical dynamic was chaotic 
and strongly diffusive. The correspondence of the observed semiclassical statistics in this chaotic 
systems and RMT means that the details of their dynamics are not relevant for the universal 
properties we have found, therefore, we expect them to hold for different underlying wave 
equations as long as they exhibit strong chaotic dynamics. 

We have found that the average conductance as a function of the chain length displays two 
regimes. First, for chains of length L < \/N, where N is the number of propagating modes in the 
leads, the dynamics is diffusive just like in a disordered wire in the metallic regime and an ohmic 
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scaling for the conductance (g(L)) = N/(L + 1) is observed. In this regime the conductance 
distribution is Gaussian and its variance grows linearly with L, but is small enough such that 
(g^ 1 ) ~ l/(<7) holds similarly to a disorder wire. Comparing the COE and CUE ensembles 
we observed a weak localization effect of magnitude similar to the usual one in the disordered 
wire. We note that the existence of a semiclassical diffusive regime is usually associated with 
disordered systems but in fact is a consequence of an underlying chaotic dynamics and can also 
be observed in periodic systems as we have shown. 

For long periodic chains, an important quantity for transport is the number of propagating 
Bloch Nb of the associated infinite system. Faure[Fau02j provided a semiclassical explicit form 
for (Nb) in diffusive systems [see equation (5.19)]; we have checked this result in a physically 
realistic waveguide and extended its validity when the unit cell possess anti-unitary symmetries 
[see equation (5.19)]. We have also linked Fame's result to the universal parametric correlation 
function [SA93a, DMS S97] IDMSS98] . Additionally, using the Machta-Zwanzig approximation 
for the diffusion coefficient, we conjectured the RMT value for (Nb) [see equation (5.36)] in 
the COE and CUE ensembles, which we checked numerically. The difference observed 5Nb = 
(Nb)coe — (Nb)cve — —0.2 shows the existence of a novel weak localization effect for the 
conductance 8g ~ SNb which extends into the asymptotic Bloch-ballistic regime [see figure [6~8] 
of a periodic chain of chaotic cavities. 

The ohmic regime breaks down at L ~ \/N and for L ^> \^N the asymptotic regime 
is reached where transport in the periodic chain is ballistic and quasi-periodic in L. This is a 
consequence of the strong resonances owing to the Bloch modes present in the associated infinite 
chain. The resonances approach the real axis as L — > oo and cluster together forming the bands 
of the infinite periodic system, thus creating Nb effectively propagative channels in the chain. 
In this regime, the transmission through these Nb channels is quasi-periodic in L. We have 
shown that the average conductance asymptotic value is (g) w 0.75(Nb). Due to the discrete 
nature of Nb the conductance distribution shows a multimodal structure in the Bloch-ballistic 



regime, with multiple peaks related to different integer realizations of Nb [see figure 6.11 . We 



showed that the conductance variance approaches an asymptotic value of order vN as L oo. 
On the other hand, the distribution of Nb resembles a Poisson distribution where the events 
correspond to the Bloch modes; however, the non-null correlation of Bloch spectra breaks down 
the exact analogy. Still, a remarkable feature is observed, namely that (Nb) /VarfiVs] approach 
a constant value as iV — > oo, which reminded us of shot noise. 

The Bloch-ballistic to diffusive transition occurring at L ~ \/N can be understood from 



several arguments. In section £12 , we noted that the distribution of the inverse decay length 
spectrum i' 1 possess a linear region in n for (Nb) < n N and for n < (Nb) the inverse 
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decay length is close to zero [see figure 6.7 due to the propagating nature of Bloch modes. The 
linear dependence of £~ = n/N inferred from the DMPK equation is the reason for the ohmic 
regime in the disordered wire. Likewise, in diffusive periodic systems, this property holds but 
with the previously named constraint which implies the ohmic regime (1/g) ~ L holds only for 
L < N/ (Nb) ~ vN. An alternative argument is that the ohmic regime should be observed 
only for chains whose diffusion time to is smaller than the unit cell Heisenberg time tn at which 
Bloch states are resolved, which leads to the same result. For L > y/N, the conductance can 
be decomposed as g = g^ + g np , with g^ the quasi-periodic asymptotic conductance which is 
a sum of the Nb Bloch modes and g np = Ae~ 2L ^ the slowest to decay non-propagating mode. 
By means of our RMT model, we have studied the pdf of t as iV — > oo and shown that under 



an appropriate scaling [see equation (6.14)] it converges to a limit distribution P(£) possessing 
an algebraic tail P(£) ~ £~ 3 for I — > oo. This allowed us to conjecture the decay (g np ) ~ L~ 2 
for the non-propagating conductance which we observed with excellent agreement in the cosine 
waveguide. 

7.1 Outlook 

There remain several open topics which are natural extensions of the work presented in this 
thesis. In this section we briefly discuss two of them that seem interesting for future research. 

7.1.1 Waveguides with Infinite Horizon 

In this work we have explicitly excluded waveguides with infinite-horizon, i.e. geometries 
allowing classical particle trajectories with free flights of infinite length. Open systems with 
strong chaotic properties possessing this property display marginally anomalous diffusion such 
that (x 2 ) = Dt log t for t oo with x t the particle position at time t and D the anomalous 



diffusion coefficient. Hence, it is obvious that Faure's result (5.19) for (Nb) must be modified 
since the normal diffusion coefficient diverge to infinity. Classically, the asymptotic billiard 
dynamics in infinite horizon systems is well understood thanks to the work of Bleher[Ble92j and 
Szazs[SV07j. They showed that the decay of correlation in the infinite horizon Lorentz Gas 
systems corresponds to the marginal case C(t) — > t^ 1 and that the random variable 

X t~ X ff7 -,\ 

Vt = ; (7.1) 
Dthgt 
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converge in distribution to a standard Gaussian law. They also gave an explicit expression for D 
in terms of the infinite horizon corridors geometry. However, the convergence of r] t is weak and 
do not imply convergence of the moments unlike the finite horizon case where even though most 
of the convergence theorems are also proved in distribution in practice it is observed that the 
finite time moments -like the finite time approximations for the diffusion coefficient- converge 



to the Gaussian moments [see discussion around equation (3.45)]. Hence, the relevance of the 
classical anomalous diffusion constant D (which can be computed analytically in some cases) 
to the semiclassical dynamics remains unclear. 

On the other hand, coming back to the derivation of Faure's expression for (Ng), there is a 



direct extension for the infinite horizon case replacing equation (5.17) appropriately. In fact, in 



the same work where (5.17) was dereived, Eckhardt [EF!C^95] gave the asymptotic dependance 
for the semiclassical variance when the classical correlation decays as C(t) — > t" 1 , namely 



login 



t 



(7.2) 



H 



Using tjj = vgjh = u^k/fiir) and conjecture (5.18), we obtain the dominant term in the 
semiclassical limit 



Hog ( 



for k — > oo. 



(7.3) 



The difference between this expression and (5.22) is numerically small in the range of energies 



we explored in this thesis. Employing a infinite horizon version of the cosine billiard, we cal- 
culated the average number of propagating Bloch modes but it was not possible to statistically 



distinguish (7.3) from (5.22) because of the large fluctuations displayed by (Ng). In order to test 



the validity of (7.3), a much deeper semiclassical region (higher energies) must be considered. 



An interesting alternative to explore would be to test (7.3) in an RMT model for periodic chains 



with infinite horizon, which could be done adding a direct process to the Circular ensembles 
of scattering matrices we have considered. This can be accomplished using the Poisson kernel, 
which is a generalization of the Circular ensembles given by the pdf 



P(5)oc|Det(l-5 t 5)|- /3 ^- 1 ^ 2 , 



(7.4) 



where 2N = i\q + N 2 with N%, N 2 the number of propagating modes on the right and left 



leads, respectively. The Poisson kernel given by (7.4) is the maximum entropy pdf satisfying 
the Circular ensembles symmetry constraints plus the condition (S p ) = (S) p = S p , where the 
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average scattering matrix S is a given sub-unit arjj^] matrix. In case S = one recovers the 
Circular ensemble with P(S) oc constant. The addition of a direct process modeling an infinite 
horizon configuration can be done defining the unit cell scattering matrix as an element of the 
Poisson kernel with 

' i 
i 



(7.5) 



where t is the transmission matrix associated to the direct process of the infinite horizon unit 
cell. 

7.1.2 Weakly disordered periodic chains 

In practice, any experimental realization of a periodic chain in the laboratory will always 
possess some degree of disorder because fabrication of trully identical unit cells is not possible. 
Thus, it would be interesting to study how weak disorder breaks down the dynamics we have 
discovered for perfect periodic chains. We have started working in this direction by considering 
an RMT model of weakly disordered time-reversal symmetric periodic chains, consisting in a 
chain of very similar but not identical cavities. We consider the case where all the cavities are 
taken from the COE ensemble; the chain ensemble is composed of cavities whose orthogonal 
scattering matrices are 

S t = U T e itH >U i = l,...,L, (7.6) 

where L G N is the chain length, U G U(N) is fixed for each realization and {Hi}f =1 is a 
realization-dependent sequence of Gaussian random real symmetric matrices. The parameter 
e is related to the disorder strength and define a particular ensemble (i.e., it does not change 
between realization of a given ensemble). The chain scattering matrix is obtained composing the 



L matrices Si by (2.13). It is clear that in the limit e = we recover the perfect periodic chain. 
Let S = U T U be an unperturbed unit cell scattering matrix. Then, all scattering matrices in 
the chain are perturbations of this original matrix, and for small e we have 

S l = S + dS i . (7.7) 



Hence, according to (4.18) and (4.19), the diference between the perturbed unit cell scattering 
matrix and the original one can be measured by fio(S — Si) = fio(dSi) ~ e . 

In order to study how disorder breaks down the transport properties of the periodic chain we 
have focused on the conductance as a function of L and e, which we denote g e {L). It is known 



* A sub-unitary matrix S is such that the eigenvalues of are < 1. 
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that in quasi-one-dimensional systems the presence of disorder of any strength will trigger 
Anderson localization, therefore for any e > we have that (g e (L)) — > as L — > oo. However, 
for sufficiently small e, we expect that in a range of L some of the periodic chain characteristics 



will still be observed. In figure 7.1, we show l/(g e (L)) for several values of e G [0,0.5] and 
N = 20. In all cases we see the ohmic linear regime for L < L max ; as e increases L max grows 
from ~ y/N as we have shown for a periodic chain to ~ N like in a disordered wire. Following 
the ohmic regime we observe different behaviors depending on the value of e. For e = we have 
that 1/ (g e (L)) — > constant for L — > oo as expected. As e > grows, we see how l/{g t (L)) starts 
diverging as L — >■ oo. We note that for small e a negative curvature region is still observed for 
L > L max before reaching exponential divergence signaling Anderson localization [see figures 



7.1 and 7.2 . This sub-exponential regime disappears for larger e suggesting that in those cases 
disorder is strong enough to erase any trace of periodicity. We have found numerically that this 
transition from weak disorder to localization occurs near e = e w 0.4 in the case considered here 
(N=20). Since fio(dSi) ~ e , we expect that as iV — > oo, this transition point moves to e — > 1. 

The previous argument shows that for weak disorder near a periodic chain there is a regime 
before localization where some reminiscent transport properties of the original system survive. 
Many interesting questions remain open, in particular about the nature of the transition from 
weak into full disorder which might be related to a percolation problem [ZS94J IIN08] . 



92 



7.1 Outlook 



Conclusions 




Figure 7.1: Inverse of the average conductance l/(g £ (L)) as a function of L for several values 
of e G [0, 0.5] and N = 20. The black line is (L + 1)/N which describes the ohmic regime and 
the dashed line is an exponential e^ L ~ 5 ^ N shown here as a guide to the eye. For e = the 
asymptotic value of 1/ (g e (L)) is a finite constant and for any e > is infinite since localization 
makes conductance null for L sufficiently large. In the latter case, it is possible to distinguish 
two different cases: for e small enough, there is still a region where l/{g e (L)) < (L + l)/iV 
with negative curvature (see figure 7.2) and for e > e this is not observed and pure exponential 
divergence after the ohmic regime is observed. We have found numerically that e ~ 0.4 for 
N = 20. 
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Figure 7.2: Derivative of the inverse average conductance as a function of L for the same e- 



realizations than figure 7.1 (the color map is the same of that figure). It can be seen that for 
small enough e there is a region with negative second derivative, i.e. where l/(g e (L)) is convex 
(negative curvature). 
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Appendix A 

Scattering matrix composition using the 

Bloch basis 



In this appendix we show an alternative method to obtain the scattering matrix of a periodic 
chain given the unit cell scattering matrix. The method consists in writing the S matrix using 
the transfer matrix Bloch spectrum. Its advantage is that it allows us to keep evanescent modes 
in the composition, hence it is more precise when the coupling between adjacent unit cells is 
not completely negligible. 

Let the eigenvectors of the unit cell transfer matrix M\ be given by F~ with associated 
eigenvalues A+, A~ where n = 1,M = . . . , N + N e , with N and N e the number of open and 
evanescent modes in the plane leads, respectively. We choose the eigenvalues A + (A - ) such that 
|A + | > 1 (|A - |<1) or 1X^1 = 1 and its associated eigenstate energy flux is right-going (left-going). 
Now, consider a periodic chain with L unit cells and let the wavefunction at the chain left end 
be 

> r and (3 = (p 



(A.l) 



n=l 



with a 



M. 



(K) L n 

T complex vectors. Therefore, applying the 



transfer matrix M = (M\) to ( A.l ), we obtain the wavefunction at the right end of the chain, 



A 2 



M 



n=l 



(A.2) 



We write the Bloch eigenvectors as 




(A.3) 



where and /„ are vectors of dimension M. We define the M x M matrices / and / with 



Scattering matrix composition using the Bloch basis 



these vectors as their columns, respectively. Then, we have that 



B 2 



r (A) L rw« 



Ri 



a 



(A.4) 



and 



r (a_)- l /- 
(A + ) L r r 

with \± the M x M diagonal matrix with the eigenvalues in its diagonal. Thus, since 



A 2 



a \ I at. 

Z3 H 2 u 



(A.5) 



A 1 
B 2 



and ip, 



out 



A 2 



(A.6) 



we can eliminate {ex., (3) from (A.4) and (A.5) to obtain 



^out = R2R1 Vin, 

which means the scattering matrix of the length L chain is given by 

/+ (a_)-*/-\/ r (A-)- L r 



(A.7) 



(A+) L f f 



(K) L r r 



(A. 
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Appendix B 

Perturbative analysis near a normal 

transfer matrix 



In this appendix we study the relation between the eigenvalues of M L and M L M L ^ when 
the transfer matrix M is close to be a normal matrix^] Let {Aj, l/\*}f =l be the spectrum of 
M and {Aj, Ar 1 }^ be the spectrum of M L M L l If M were normal we could map these spec- 
trums such that Aj = |Aj| 2L , hence the transfer eigenvalues Tj = 4 (2 + Aj + Aj" 1 )^ 1 associated 
to propagating Bloch modes would be not dependent on L and equal to one, so for long chains 
g = Nb- Generically, the transfer matrix is not normal, leading to a more involved relation 



between the sets of Aj and Aj as a function of L [see discussion around (2.35)]. 



We consider a general unit cell transfer matrix M constrained only by (2.19), which implies 
S-pseudo-unitarity and no further symmetries. Without loss of generality we can always write 
the L-unit-cells transfer matrix in the form 

M L = Pe iDL p-\ (B.l) 

where the matrix D is diagonal and composed of N complex conjugate pairs {#j, 9*}f =l . Of 
course, e lDL is the diagonal matrix with M eigenvalues Aj = e l9lL . The 2A?# propagating Bloch 
modes eigenvalues (|Aj| = 1) are given by real D elements. The matrix M is normal if and only 
if its eigenvectors matrix P is unitary. We start the perturbative expansion assuming that 

P = P + tP 1 (B.2) 

with Pq unitary and e C 1, i.e. M is close to be normal. We now proceed to expand 

M L M Lt = p e iDL p -lp-l] e -iD*Lp] (B.3) 



*A matrix is said to be normal if it commutes with its Hermitic conjugate or equivalently if its eigenvectors 
form an orthogonal set. 



Perturbative analysis near a normal transfer matrix 



to order e. First, we take the inverse of P, 

P 1 = (P + eP,)- 1 = Pj - ePjPiPj + 0(e 2 ) (B.4) 

Then, 

p-ip-it = 1 _ e ^ptp Q + ptp^ + (g2) (B.5) 
where we have used that P Pq = 1. Finally, we obtain 

M L M Lt = P [H + eJETi] P f + £>(e 2 ) (B.6) 

with H = e *(«-D*)£ = e-'W and 

H = ifoP/Po + PjiWo - e ,;DL (P/Po + PjPi) e" iD * L . (B.7) 

Note that H = Hq + cHl is Hermitic so the eigenvalues of M L M Lj[ are equal to the eigenvalues 
of H. Hence, in order to obtain this spectrum we can use the well known perturbation theory 
for time-independent Hamiltonians[Sak94j. We want to calculate the spectrum given by the 
eigenvalue problem 

Hcj) = A<f>. (B.8) 

As usual, we set A = A + eA 1 . Since our unperturbed Hamiltonian H is already diagonal, 
we have that A° = e ~ 2:,m(ei)L , for i — 1, . . . ,2N, with 9 %i = Q* 2i _ x . Note that the unperturbed 
spectrum {A°, (A°)~ 1 }^ 1 is degenerate with the 2Nb eigenvalues related to propagating Bloch 
modes equal to one. Also, it is clear that the transmission eigenvalues Tj associated to non- 
propagating modes decays exponentially so we are interested in calculating only the correction 
to the degenerate set of Bloch transmission eigenvalues. 

The first order correction A 1 to the 2A#-degenerate eigenvalues A? = 1 is given by the 
solutions of the eigenvalue problem 

h[ 2Nb) ^ = A] fa (B.9) 

with JJ^ 2Nb ^ the projection of Hi to the degenerate eigenspace. If we choose the 2Nb unit 
eigenvalues at i — 1,...,2N B then h[ 2Nb ^ is the upper left 2N B x 2N B block of Hi. The 
transmission eigenvalues associated to the degenerate Bloch subspace are given to dominant 
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order in e by 



(B.10) 



We can illustrate this result by means of the simple case Nb = 1. When there is only one 



propagating Bloch mode, it is trivial to obtain the solutions to (B.9), which gives 



A± = l±2\v[v 2 sin [(e 1 -9 2 )L]\, 



(B.ll) 



where P = [v\ v 2 ] with v\v 2 ~ 0{t). In this case, the transmission eigenvalue is 



T x = l-\v\v 2 sin[(^-0 2 )L]| 2 . 



(B.12) 



This shows how the non-orthogonality of M eigenvectors breaks A degeneracy creating a peri- 
odic dependance of T\ on the chain length L. 
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Appendix C 

Analytical calculation of (Ng) for N = 1 



In this appendix we calculate analitically the average number of propagating Bloch modes 
(Nb) for a COE and CUE periodic chain model with one open mode in the plane leads, this 
is N — 1. In this case the scattering and transmission matrices are 2x2 allowing an explicit 
analitycal solution of (Nb)- 



From definition (2.18), of the transfer matrix M as a function of the scattering matrix S, 
we have that for N = 1, 

' t - rr'/t 1 r'/t 1 
-r/t' 1/f 



M 



(C.l) 



where t, r (t',r') are the left-to-right (right-to-left) transmission and reflection coefficients, re- 



spectively. They are complex numbers and define the scattering matrix [see section 2.1.1 




(C.2) 



As we have shown in section 2.1.2, the number of propagating Bloch modes Nb in a periodic 



chain with unit cell S-matrix given by (C.2) is equal to the number of iVf-matrix eigenvalues A± 



with unit modulus. In the present case, we can solve the M eigenvalue problem algebraically 
to obtain 

rr' + tt r 



A 



± 



1 



2f 




(C.3) 



Since there is only one open channel in the leads, the Bloch spectrum consists of a pair of 
states corresponding to modes traveling in both chain directions. As we will see explicitly later, 
lA-l -1 = IA+I. Therefore, we have that 



(N B (N = 1)) = P[A± G S 1 ], 



(C.4) 



where P[\± G S 1 } is the probability that |A±| = 1 under the appropriate RMT measure. We 



start with the CUE periodic chain ensemble. As discussed in section 4.3, the ensemble of 



Analytical calculation of (N b ) for N = 1 



scattering matrices S taken from the CUE ensemble are distributed according to the Haar 
measure of the unitary group U(2iV). In general, any matrix of U(2iV) can be represented using 
Euler angles as a composition of rotations in U(2) [ZK94J. For N = 1, a general S matrix is 
described by four independent parameters as 



>1,<P2,<P3 



-101 



(C.5) 



with £ G [0, 1] and <pi, <f) 2 , 4>3 G [0, 2ir). The Haar measure in this parametrization is simply 

1 



fJLU(2)(dS) 



(27T)' 



)i a<p 2 acpz- 



Now, in order to compute P[X± G S 1 ], we rewrite (C.3) using (C.5) as 



A ± = e <(*/2-fc) (a ± V-l + a 2 ) 



with 



Thus, since a G 



a 



1-rr' + tt' 



sin 03 



I A 



± 



1 if and only if |a| < 1. Hence, we obtain 



P[\± eS 1 ]= [ 

Vis 



fiu(S){d8) 

I sin (fs\ 2 <l-( 

2tt /•!-! sin</> 3 | 2 



(2 



7T 



^0 



The integral in (C.10) can be easily computed and by (C.4) leads to 



(iV fl (JV=l)> 



1 



CUE 



(C.6) 



(C.7) 



(C. 



(C.9) 
(C.10) 



(C.ll) 



which is an exact result. This means that in a CUE periodic chain with one open channel the 
events of a Bloch state being propagative or evanescent are equiprobable (and complementary). 



We note that the asymptotic result obtained in section 5.2, namely (N b (N)) C ve = 

y/N/^ for 



N — > oo, in this case gives y 1/7T = 0.56. 

We now turn to the COE periodic chain. As we have shown in section 4^3, an S matrix 
in this ensemble can be constructed by taking U G U(2N) distributed according to the Haar 
measure /iu(2N) an d setting S = U T U. Since for N = 1 we have an explicit parametrization for 
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Analytical calculation of (N b ) for N = 1 



U given by (C.5), we can write a general symmetric unitary matrix in COE as 



1 -2i v '{(l-«)siii(^ 1 +^) «e-**' + (l-fl<** ' 



where again £ G [0, 1] and 01,02,03 € [0, 27r) are distributed according to the uniform measure 



(C.6). In this case, the eigenvalues of M are 

A± = a' ± V-l + a' 2 (C.13) 

with 

^= o Si ^^ csc(0 1 + 2 ) (C.14) 

a real function of the parameters. Therefore, we have again that |A±| = 1 if and only if \a'\ < 1, 
which implies 

P[X± eS x ]= [ VumidS) (C.15) 

J|sin203| 2 <2 A /?(l-OI sin (<£i+<fe) 

= -^r [ [ arcsin (2Ml - £) I sin flf) d0 . (C.16) 
2tt j Jo Jo ^ J 

This integral can be evaluated numerically to obtain 

(N B {N = 1))coe = 0.3634. (C.17) 
This result is quite close to the value given by the asymptotic expression [see section |5.2| 



(N B (N)) COE = V N / n ~ °- 2 for N -> °°> wnich in this 

case gives ■s/l/^- 0.2 = 0.3642. 
The present procedure can in principle be generalized to arbitrary N, however for larger 
dimensions the method becomes increasingly involved and impractical, since eigenvalue char- 
acteristic polynomials are not solvable analytically and, in addition, the direct evaluation of 
integral expressions is not feasible in high dimensional spaces. 
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